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CHAPTER 1 



Introduction and overview 



1.1 A brief overview on magnetic fields in the cosmos, their 
observation, their origin and particle physics 

Magnetic fields are most frequently found in various length scales and various positions 
of our universe. The average magnitude of the Earth's magnetic field is about lO^^G. This 
magnetic field encircles the Earth and produces a magnetosphere surrounding it. The field 
is not uniform. Charged elementary particles from the cosmos sometimes gets trapped in 
this magnetic field and hover over our Earth in a belt shaped region which is popularly 
called the Van Allen belt. Also the solar storms produce electromagnetic disturbances on 
the magnetic field of the Earth and the particles trapped in the Van Allen belt, resulting 
in 'aurora borealis', one of the strangest things seen by human beings. In the solar black 
spots magnetic fields can go up to lO^G. Very high magnetic fields are found in the cores 
of supernovas, and this magnetic field remains frozen on the proto-neutron star remnant 
whose surface magnetic field has been measured to be lO^^G. Not only the planets and the 
stars have magnetic fields associated with them there are magnetic fields in the galactic 
and intergalactic mediums. In Milky Way the average galactic magnetic field is of the order 
of 3 — 4G with a length scale of a few Kilopersecs. Average intergalactic magnetic fields 
can be of order 10~^^G P with a length scale of about 1 Megaparsec. Except gravitational 
interactions which becomes important at astronomical distances (some Megaparsecs) at 
average energies (much lesser than the Planck scale lO^^Gev) the only other field which is 
omnipresent in the universe in various shapes, sizes and magnitudes is the magnetic field. 

The main observational tracers of galactic and extra-galactic magnetic fields are jU |2I 
the Zeeman splitting of spectral lines, the intensity and polarization of synchrotron emis- 
sion from free relativistic electrons and Faraday rotations of polarized electromagnetic 
radiation passing through an ionized medium. The Zeeman splitting is too small to be 
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useful outside our own galaxy. The synchrotron emission method and Faraday rotation 
measurements require an independent estimate of the local electron density. If the mag- 
netic field to be measured is far away one relies on Faraday rotation. The agreement 
generally found between the strength of the field determined by Faraday rotation tech- 
nique and that inferred from the analysis of the synchrotron emission in relatively close 
objects gives reasonable confidence on the reliability of the first method also for far away 
systems. 

The origin of the magnetic fields observed in the galaxies and in the cluster of galaxies 
is unknown. This is an outstanding problem in modern astrophysics and historically it was 
the first motivation to look for a primordial origin of the magnetic fields jS] . The general 
trend is to use magnetohydrodynamic methods to amplify very weak magnetic fields into 
the /iG fields generally observed in galaxies. Today the efficiency of such a mechanism for 
production of magnetic fields is in question from new observations of magnetic fields in 
high redshifted galaxies. Furthermore, even if the magnetohydrodynamical calculations 
are taken into account, the origin of the seed fields which initiated the process has still to be 
identified. It is understood that somewhere elementary particle physics will enter the scene 
to explain the formation of the elusive seed fields. With the cosmological observational 
evidences, as cosmic micro wave background radiation (CMBR) and nucleosynthesis data, 
big bang theory is a cosmological reality now and so is standard big bang cosmology. In 
near future perhaps we will have a standard mechanism for production of the seed fields 
coming out from big bang cosmology. 

In the previous paragraph we discussed how models of particle physics can be used 
to produce the seed fields. But the magnetic fields thus produced can themselves af- 
fect various elementary particle interactions. Magnetic fields can enhance the scattering 
cross-sections or decay rates from their vacuum values. The action of magnetic fields on 
elementary particle interactions can produce non trivial results and can explain some of 
the interesting observations in astrophysics. The main portion of this thesis deals with 
elementary particle interactions, of which the neutrino is of prime concern, in presence 
of a uniform background magnetic field. To put the theme in perspective the following 
section gives an overview of the nature of work done in this field by previous workers and 
the main points around which this thesis evolves. 
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1.2 Some well known observations, results and motivation for 
this thesis 



Neutrinos have no electric charge. So they do not have any direct couphng to photons 
in any renormahzable quantum field theory. The standard Dirac contribution to the 
magnetic moment, which comes from the vector coupling of a fermion to the photon, is 
therefore absent for the neutrino. In the standard model of electroweak interactions, the 
neutrinos cannot have any anomalous magnetic moment either. The reason is: anomalous 
magnetic moment comes from chirality-fiipping interactions ipa^^yipF^'^ , and neutrinos 
cannot have such interactions because there are no right-chiral neutrinos in the standard 
model. Other particles as the electron and muon have normal and anomalous magnetic 
moments and so can interact with external magnetic fields. 

Inclusion of neutrino mass naturally takes us beyond the standard model, where the 
issue of neutrino interactions with a magnetic field must be reassessed. If the massive 
neutrino turns out to be a Dirac fermion, its right-chiral projection must be included in 
the fermion content of the theory, and in that case an anomalous magnetic moment of a 
neutrino automatically emerges when quantum corrections are taken into account. In the 
simplest extension of the standard model including right-chiral neutrinos, the magnetic 
moment arises from the diagrams in Fig. 11.11 and is given by j3] 



where is the mass of the neutrino and /i^ is the Bohr magneton. If, on the other hand, 
neutrinos have Majorana masses,^ i.e., they are their own antiparticles, they cannot have 
any magnetic moment at all, because CPT symmetry implies that the magnetic moments 
of a particle and its antiparticle should be equal and opposite. However, even in this case 
there can be transition magnetic moments, which are co-efficients of effective operators 
of the form ip-^^afj_uip2F'^'' , where ipi and 1IJ2 denote two different fermion fields. These will 
also indicate some sort of interaction with the magnetic field, associated with a change of 
the fermion fiavour. The question of the neutrino magnetic moment assumed importance 
when it was suggested that it can be a potential solution for the solar neutrino puzzle 
[HI 13 El- A viable solution required a neutrino magnetic moment around 10^^^ ^b, orders 
of magnitude larger than that given by Eq. (|l.ip , knowing that the neutrino masses cannot 

detailed discussion on Dirac and Majorana masses of neutrinos is present in [S]. 
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Figure 1.1: One-loop diagrams that give rise to neutrino magnetic moment in standard model 
aided with right-handed neutrinos. The lines marked v are generic neutrino lines, whereas 
those marked t are generic charged leptons. The external vector boson line is the photon. In 
renormalizable gauges, there are extra diagrams where any of the W lines can be replaced by 
the corresponding unphysical Higgs scalar. 

be very large. At present the situation has changed and neutrino magnetic moment is not 
a lucrative object to look for to explain the solar neutrino data. The Sudbury neutrino 
observatory results [9, indicate no deficiency of left handed neutrinos in neutral-current 
interactions. 

If we leave out non-standard neutrinos and take its magnetic moment to be zero, 
then also magnetic fields can alter their properties. This is the prime concern of the 
present thesis. One of the simplest things that gets affected is the neutrino self-energy 
in a magnetic field. The magnetic field enters the scene through interaction with the 
charged leptons in the loops of the self-energy diagram. Erdas and Feldman calculated 
the dispersion relation of neutrinos in a uniform background magnetic field. To lowest 
order in the magnetic field B they obtained the dispersion relation |10j 



Here f2 and q are the neutrino energy and 3-momentum. is the magnitude of that 
component of the neutrino 3-momentum which is perpendicular to the magnetic field 3- 
vector. m is the mass of the charged fermion which has the same flavour as that of the 
neutrino, My^ the W boson mass and g the weak coupling constant. The origin of such a 
term like will be discussed in chapter El The dispersion relation reduces to the usual 
dispersion relation of a zero mass particle when we put i3 = 0. For strong magnetic fields. 
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the dispersion relation has been calculated more recently by Elizalde, Ferrer and de la 
Incera 

There are other processes as neutrino neucleon scattering which are also affected by 
magnetic fields. The charged current interaction Lagrangian involving neutrinos and 
nucleons is given by 

=^nt = ^Gfi ['0(^)7'' L?/^(^,)] [V'(p)7m(G'v + G'A75)V'(n)] , (1-3) 

where L = |(1 — 75) and = Gpco^Qc-, Oc being the Cabibbo angle, and GaIGv = 
—1.26. This can be used to find the cross-section for various neutrino-nucleon scattering 
processes. 

The above interaction Lagrangian is most vital for processes involving neutrinos or 
antineutrinos. In a star, when such reactions occur, the final neutrino or the antineutrino 
escapes and the star loses energy. Such processes are collectively known as URCA pro- 
cesses, named after a casino in Rio de Janeiro where customers lose money little by little 
[12]. These processes are, 

n — > p + e~ + z7 , (1.4) 
p + e~ — > n + 1/ , (1.5) 
n + ^ p + V . (1.6) 

The calculations of the cross-sections and decay rates in presence of a magnetic field uses 
a benchmark value of the background field: 

B = — 

e 

= 4.41x10^=^0, (1.7) 

where m is the electron mass, and e is the charge of the proton. is sometimes called 
the 'critical field' although nothing that critical happens at this field strength. The only 
condition which defines the critical field is that when the magnitude of the magnetic field 
reaches the value Be the electron cyclotron frequency equals its rest mass in the natural 
units. If the magnitude of the magnetic field is above the critical field the electronic 
wave functions get considerably modified. It is assumed henceforth that while discussing 
scattering or decay processes the magnitude of the external fields will be around Be- If 
the magnitude of the field is comparable to m^/e then the proton wave functions can also 
get affected. The calculations in this thesis will not involve such high fields. 



The rate of the neutron beta decay process in a magnetic field was calculated by various 
authors. An early paper by Fassio-Canuto |13i derived the rate in a background of degen- 
erate electrons. The work considers exact wave functions of the electrons in a background 
of a uniform magnetic field. Contemporary papers by Matese and O'Connell ^1 
derived the rate where the background did not contain any matter, but included the 
effects of the polarization of neutrons due to the magnetic field. Protons and neutrons 
were assumed to be non-relativistic in the calculations. Further, the magnetic field was 
assumed to be much smaller than m^/e so that its effect on the proton wave function 
could be neglected. Various calculations of the other processes exist in the literature. 
Some calculations take the background matter density into account [121 HZl UHl CHi 120] , 
some include magnetic effects on the proton wavefunction as well j^U]- In all of these 
calculations the cross-section or the decay rate is sensitive to the direction between the 
neutrino 3-momentum and the external uniform magnetic field. 

Neutrinos can also interact with photons, but this interaction is an effective one me- 
diated by virtual charged particles. In a magnetic field the propagators of these charged 
particles are affected by the external field. Previously it has been shown how the neu- 
trino self-energy is modified in an external magnetic field. The properties of photons also 
get modified in an external uniform magnetic field due to their interaction with virtual 
charged particles. It was shown by Adler |^ that when the background magnetic field 
magnitude is around many interesting things can happen. The photon dispersion re- 
lation in presence of a magnetic field changes from that in vacuum |A;| = uj. There are 
actually two different indices of refraction corresponding to the two photon propagation 
modes. The modes are linearly polarized, with the magnetic field of the mode either 
parallel or perpendicular to the plane containing the external magnetic field and the 
direction of propagation. Most of the calculations done by Adler uses the effective Euler- 
Heisenberg Lagrangian |j22j. Later Schwinger, while trying to derive the Euler-Heisenberg 
Lagrangian, derived the form of the propagator of the charged particles in an uniform 
background magnetic field p3]. After Schwinger's work most of the calculations, as the 
self-energy of the photon in a magnetic field, are done using the Schwinger propagator. 

Two second rank tensors play an important part in building up the effective neutrino 
photon vertex function. They are the photon self-energy 11^,^ and another tensor 11^^. To 
one loop the tensors 11^^, and Il^i, are similar except one point. In Tl^i, both the vertices, 
corresponding to the two tensor indices, are of vector type whereas in n|^^ one of the 
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vertices is of the axial- vector type |21I. 11^^^ is called the axialvector- vector amplitude in 
this thesis. In a magnetic field background 11^,^ was first calculated by Tsai [23]. In a 
different work DeRaad, Milton and Hari Dass [20] calculated 11^^,. Both these calculations 
were to one loop but to all orders in the external magnetic field. Using the expression of 
the above two tensors in a magnetic field the rate of the neutrino Cherenkov process, 



where g, g' and k are the four momenta of the particles as shown, was calculated by Raffelt 
and loannisian |2Zl I2H1- The neutrino Cherenkov process is kinematically forbidden in 
vacuum. But in a magnetic field as the photon dispersion relation changes so this process 
becomes viable. The photon self-energy U^^, [23 1^0] and 11^^^ jSH E21 IHH] both have been 
calculated in a magnetized medium to one loop. 

In this thesis the main emphasis has been given on neutrino scattering processes in 
presence of a magnetic field and the electromagnetic vertex of neutrinos in a magnetized 
medium. 11^^, which is related to the electromagnetic vertex of neutrinos, has been 
calculated in the background of a magnetized medium. Using the expression of 11^^ 
attempt has been made to find out a specific limit of this quantity, which gives the 
'effective electric charge' of the neutrino. A general form factor analysis of 11^^ for various 
backgrounds, as in vacuum, thermal medium, magnetic field and a magnetized medium 
is also given in chapter |3] Neutrino interactions in presence of a classical background 
magnetic field is a subject which has been studied for a long time and most of the topics 
mentioned above are briefiy discussed in the review by Bhattacharya and Pal 

1.3 Notations and conventions 

Before going into the second chapter, this section summarizes the mathematical no- 
tations and conventions we use throughout this thesis. In future the magnitude of the 
uniform background magnetic field will always be denoted by the symbol B and the mag- 
netic field vector pointing in the z direction will be given by 



where z is the unit vector along the z-axis. 

The first thing to note is that a uniform classical magnetic field chooses a certain pre- 
ferred direction in space. Consequently the Lorentz invariance of the system is restricted, 
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B = zB, 




arbitrary boosts will not preserve a pure magnetic field. A boost along the magnetic field 
direction or a rotation about the external field direction, are the only Lorentz transfor- 
mations allowed now. As a result the 4- vector structure breaks down into a perpendicular 
part and a parallel part. If a'^ is a 4-vector, then 

af = (a°,0,0,a=^), (1.10) 
= (0, a\ a^, 0) , (1.11) 

such that 



Also in our convention 



a^^^ajf + a^^. (1.12) 



9^u = glu + gf^u^ (1-13) 



where 



(1,0,0,-1), (1.14) 

= (0,-1, -1,0). (1.15) 

Some times instead of ay and aj^ use has been made of a'^n and a^^ , they are the same 
things written in two different ways. 

Frequently we will come across terms as a| and a\. These terms stands for 

aJ^iaV, (1.16) 

= (aT-(a^)^ (1.17) 

ai = (1-18) 

= {a^f + {a^)\ (1.19) 

such that 

^a\-a\. (1.20) 

Sometimes the equations will involve terms as (a • &)||, where a and h are two 4- vectors. 
This dot product stands for 

(a • 6)|| = (a°6^ - a^6°) = -(0063 - as&o) , (1-21) 
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which imphes 

(a ■ = -{b ■ 5)j| , 

and result 

(a ■ a)|| = . 

1.4 Scheme of things to come 

This last section is about the scheme of the things to come. Chapter |21 focuses on the 
quantum field theoretical description of scattering processes which involves charged parti- 
cles, in presence of a background uniform magnetic field. The field theoretical techniques 
are then applied to find out the cross-section of the inverse beta decay process. Using the 
expression of the inverse beta decay cross-section an attempt has been made to explain 
the puzzle regarding the high velocities of pulsars. 

Chapter 01 discusses about the charged particle propagators in an uniform background 
magnetic field using Schwinger's method ''^\. The Schwinger propagator is introduced 
and then by using the methods of quantum statistical field theory it is expressed in a 
thermal background. A general analysis on neutrino self-energy in a magnetic field and 
in a magnetized medium follows. The chapter ends with a discussion on the effects of a 
background magnetic field on neutrino oscillation phenomenology. 

Chapter m deals with the topic of neutrino photon interactions. It is shown that the 
photon vacuum polarization tensor U^i, and the axialvector-vector amplitude 11^^ are 
essential ingredients of the neutrino electromagnetic vertex. Some comments on U^,^ and 
a discussion on 11^^, is presented in this chapter. The chapter ends with an expression of 
the effective electric charge of neutrinos in a magnetized medium. 

The last chapter concludes by summarizing the points which are discussed in the 
various chapters of this thesis. The thesis ends with a set of appendices containing some 
detailed calculations which were put at the end so as not to block the readers mind with 
lengthy calculations. 
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CHAPTER 2 



Inverse beta decay in a constant 
background magnetic field 

2.1 Introduction 

The interactions of elementary particles show novel features when they occur in non- 
trivial backgrounds. Study of particle propagation in matter has proved pivotal in the 
understanding of the solar neutrino problem. A typical example of this is the neutrino- 
electron scattering in presence of a medium which culminated in the now famous MSW |SSl 
EE] effect. The cross-section for the elastic neutrino-electron scattering in presence of a 
magnetic field has been studied j^Z]. Bezchastnov and Haensel |3B.^ has also calculated 
the neutrino-electron scattering cross-section in presence of a magnetized medium. 

There are processes related to the neutrino-electron scattering, obtained by crossing, 
which contain a neutrino-antineutrino pair in the final state. For example, one can have 
the pair annihilation of electron and positron into neutrino-antineutrino pair, or the neu- 
trino synchrotron radiation: 

e"^e^ + z/ + z7. (2.1) 

This is similar to a the normal synchrotron radiation reaction, the difference being that a 
neutrino-antineutrino pair is produced instead of a photon. It should be noted that this 
process cannot occur in the vacuum. However, in the presence of a background magnetic 
field and background matter, the dispersion relation of the electron changes so that it 
becomes kinematically feasible. These processes provide important mechanism for stellar 
energy loss, and the rates of these processes have been calculated jSHl HOI E] . The reverse 
of these processes are important for neutrino absorption, and have also been studied • 
Moreover studies of scattering cross-sections and decay rates in background magnetic 
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fields are important since stellar objects like neutron stars are expected to possess very 
high magnetic fields, of the order of 10^^ G or higher. Analysis of these processes might 
be crucial for obtaining a proper understanding of the properties of these stars. As 
the preferred direction of the background magnetic fields break the isotropy of space the 
scattering cross-sections of various processes calculated in such a background also becomes 
anisotropic. This anisotropy of the cross-sections of processes involving neutrinos as initial 
particles or final particles can give rise to an overall asymmetric emission of neutrinos from 
a newly born neutron star and as a result the star may get an overall momentum thrust. 

As discussed in chapter ^ there are various possible ways to introduce the external 
uniform magnetic field in the calculation of processes including elementary particles. In 
this chapter we explore the exact solutions of the Dirac equation in the presence of such 
an external field, and then apply the results to obtain the physically interesting numbers 
as scattering cross-sections IH7] . The initial sections of this chapter develops a consis- 
tent quantum field theoretical method for calculations of cross-sections in a background 
magnetic field. The discussions up to section IT^ is general, the results discussed can be 
appropriately used in any scattering cross-section calculation in a background magnetic 
field. The particularities of the inverse beta decay process cross-section calculation enters 
from section 12.41 Based upon the form of the scattering cross-section, section 12.51 eluci- 
dates the point about asymmetric emission of neutrinos from a proto-neutron star. The 
chapter concludes by summarizing the ideas exposed. 

2.2 Solutions of the Dirac equation in a uniform magnetic field 

For a particle of mass m and charge eQ, the Dirac equation in presence of a magnetic 
field is given by 



where a and /? are the Dirac matrices, and A is the vector potential. In our convention, 
e is the positive unit of charge, taken as usual to be equal to the proton charge. 
For stationary states, we can write 






(2.3) 
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where (j) and x ^^re 2-component objects. We use the Pauh-Dirac representation of the 
Dirac matrices, in which 



cr 
a- 



1 
1 



(2.4) 



where each block represents a 2 x 2 matrix, and cr are the Pauh matrices. With this 
notation, we can write Eq. ()2.2|) as 



{E — ■m)(f) 
{E + m)x 



Ehminating x, we obtain 



(7-{-tV -eQA)x. 
cr ■ {-iV - eQA)(f). 



a ■ {-iV - eQA) 



(2.5) 
(2.6) 



(2.7) 



We will work with a constant magnetic field B. Without loss of generality, it can be taken 
along the ^-direction. The vector potential can be chosen in many equivalent ways. We 
take 



Ao = Ay = A, = 



A^ 



-yB. 



With this choice, Eq. ()2.7|) reduces to the form 



d 



V + {eQB)y - eQB{2ty— + a-, 



(2.8) 



(2.9) 



Here is the diagonal Pauli matrix. Noticing that the co-ordinates x and z do not 
appear in the equation except through the derivatives, we can write the solutions as 



(2.10) 



where f{y) is a 2-component matrix which depends only on the y-coordinate, and possibly 
some momentum components, as we will see shortly. We have also introduced the notation 
X for the spatial co-ordinates (in order to distinguish it from x, which is one of the 
components of X), and X^ for the vector X with its y-component set equal to zero. In 
other words, p ■ X^ = p^x +PzZ, where and Pz denote the eigenvalues of momentum in 
the X and z directions.^ 



^It is to be understood that whenever we write the spatial component of any vector with a lettered 
subscript, it would imply the corresponding contravariant component of the relevant 4-vector. 
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There will be two independent solutions for f{y), which can be taken, without any 
loss of generality, to be the eigenstates of with eigenvalues s = ±1. This means that 
we choose the two independent solutions in the form 

Since a^fs = sfs, the differential equations satisfied by Fg is 

'^'^^ (eQBy + p^fF, + {E^ -m^-pl + eQBs)F, = , (2. 12) 



which is obtained from Eq. ()2.9|) . The solution is obtained by using the dimensionless 
variable 

Px 



f = V^MSh/+^ . (2.13) 



which transforms Eq. (j2.12j) to the form 

rf2 



0, (2.14) 



where 

E'^-m'^-pl + eQBs 



(2.15) 



e\Q\B 

This is a special form of Hermite's equation, and the solutions exist provided = 2z/ + 1 
for z/ = 0, 1, 2, ■ ■ -. This provides the energy eigenvalues 

E'^ = m'^+pl + (2z/ + l)e\Q\B - eQBs , (2.16) 

and the solutions for Fg are 

N^e^^'/'H^iO^m, (2.17) 

where Hi, are Hermite polynomials of order z/, and N,y are normalizations which we take 
to be 

1/2 



/e\Q\B 



We stress that the choice of normalization can be arbitrarily made, as will be clarified 
later. With our choice, the functions I^, satisfy the completeness relation 



J2 UOU^^) = ^e\Q\B - e.) = 5iy - y.) , (2.19) 

V 
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where is obtained by replacing y by y^, in Eq. ()2.13|) . 

So far, Q was arbitrary. We now specialize to the case of electrons, for which Q = —1. 
The solutions are then conveniently classified by the energy eigenvalues 

El = m^ +pl + 2neB, (2.20) 

which is the relativistic form of Landau energy levels. The solutions are two fold degen- 
erate in general: for s = l, z/ = n — 1 and for s = —1, v = n. In the case of n = 0, 
only the second solution is available since v cannot be negative. The solutions can have 
positive or negative energies. We will denote the positive square root of the right side by 
En- Representing the solution corresponding to this n-th Landau level by a superscript 
n, we can then write for the positive energy solutions. 



■4-1(0' 



(2.21) 



For 77, = 0, the solution /+ does not exist. We will consistently incorporate this fact by 
defining 

/-i(z/) = 0, (2.22) 

in addition to the definition of /„ in Eq. ()2.17|) for non-negative integers n. 

The solutions in Eq. ()2.21|) determine the upper components of the spinors through 
Eq. ()2.1Up . The lower components, denoted by x earlier, can be solved using Eq. ()2.6p . 
and finally the positive energy solutions of the Dirac equation can be written as 

e-^P-^^?7,(2/,n,p^), 



(2.23) 



where denotes the space-time coordinate. And Us are given by 



U+{y,n, 







Pz 



En + m 



V2neB 



V En + 



m 



4(0 



U-iy,n,: 



V 





4(0 

V2neB , 
En + m 

Pz 

En + m 



(0 
4(0 



. (2.24) 



A similar procedure can be adopted for negative energy spinors which have energy 
eigenvalues E = —En- In this case, it is easier to start with the two lower components 
first and then find the upper components from Eq. (j2.5|) . The solutions are 

e*P-^^V;(i/,n,p^), (2.25) 
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where 



/ Pz 



En + m 



(0 



V 



V2neB 
En + m 





Pz 



J 



En + m 




\ 



(2.26) 



where ^ is obtained from by changing the sign of the pa;-term. 

For future use, we note down a few identities involving the spinors which can be 
obtained by direct substitutions of the solutions obtained above. The details of the 
calculation has been worked out in appendix |X1 The spin sum for the [/-spinors is 



Y,Us{y,n,j^)Us{y^,n,p^) 

s 

1 



2{En + m) 



X 



|m(l + a^) +^ii -|^||75|/„_i(0/„_i(^^ 

+ {m(l - a,) + + ^1,75} /„(04(e.) 
-72^(71 - z72)/n(04-i(e.) 
-72^^(71 + z72)/n-i(04(e. 



(2.27) 



where the above notations for the 4- vectors and the || and _L notations has been defined 
in chapter n and = ^7^7^. The symbol = p'^73 +p^7o- Similarly, the spin sum for 
the l^-spinors can also be calculated, and we obtain 



Pv{y,y*,n,j\) 



Y,Vs{y,n,i\,)Vs{y,n,i^) 

s 

1 

2{En + m] 



X 



{-m(l + (12) +1^11 - 1^1175} /„_i(^)/„_i(^* 

+ {-"^(1 - (^z) + h t?||75} 4(04(1*) 
+72^^(71 - «72)4(0/n-i(e.) 



+V2^(7i + «72)4_i(e)/„(e*) 



(2.28) 
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2.3 The fermion field operator 



Since we have found the solutions to the Dirac equation, we can now use them to 
construct the fermion field operator in the second quantized version. For this, we write 



m) = E E 



s=± n=0 



D 



f,{n,p^)e-^P-''^U,{y,n,p^) + fj{n,p,)e^^-''^Vs{y,n,p^) (2.29) 



Here, /s(n, p^) is the annihilation operator for the fermion, and fliji^p^) is the creation 
operator for the antifermion in the n-th Landau level with given values of Px and Pz. The 
creation and annihilation operators satisfy the anticommutation relations 



(2.30) 



and a similar one with the operators / and all other anticommutators being zero. The 
quantity D appearing in Eq. ()2.29|) depends on the normalization of the spinor solutions, 
and this is why the normalization of the spinors could have been chosen arbitrarily, as 
remarked after Eq. ()2.18|) . Once we have chosen the spinor normalization, the factor D 
appearing in Eq. ()2.29|) is however fixed, and it can be determined from the equal time 
anticommutation relation 



6'(x - x^ . 



(2.31] 



Plugging in the expression given in Eq. ()2.29|) to the left side of this equation and using 
the anticommutation relations of Eq. ()2.30|) . we obtain 

'^P{X),^P\X.)] = EE / (e-^^^(^-^*)e-^^^(-*)f/.(y,n,p^)f/t(y.,n,p^) 

s n •' ^ 

+e^"^^^-^*)e^^^(^-^*)K(2/,ri,p^)K^(|/.,n,F^)). 

(2.32) 

Changing the signs of the dummy integration variables px and pz in the second term, we 
can rewrite it as 

V^(X),V^t(x,)l^ = ^^ e-^^^(—*)e-^^^(— *)([/.(!/, n,p^)^t(^^,^,p^) 



(2.33) 
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Using now the solutions for the U and the V spinors from Eqs. ()2.24|) and ()2.26|) . it is 
straight forward to verify that 



J2 [Us{y, n, p^)Ul{y^, n, p^) + Vs{y, n, -p^)Vj{y^, n, -p^ 

pl + 2neB\ 
(E„ + m)2j 



r)2 _|_ 2neB\ r 



(2.34) 



where 'diag' indicates a diagonal matrix with the specified entries, and C, and .^^ involve 
the same value of px- At this stage, we can perform the sum over n in Eq. ()2.33|) using the 
completeness relation of Eq. ()2.19p . which gives the 5-function of the y-coordinate that 
should appear in the anticommutator. Finally, performing the integrations over and 
Pz, we can recover the 5-functions for the other two coordinates as well, provided 

using the expression for the energy eigenvalues from Eq. (|2.2Up to rewrite the prefactor 
appearing on the right side of Eq. (j2.34j) . Putting the solution for D, we can rewrite Eq. 
as 



oo „ 



dpxdpz En + m 



s=± n=0 • 
X 



27T V 2E„ 

[fs{n, p^)e-^^-^^f/,(i/, n, p,) + fUn, p,)e'^-''W,{y, n, p,)] . (2.36) 

The one-fermion states are defined as 

\n,p,) = Cpin,p,)\0) . (2.37) 

The normalization constant C is determined by the condition that the one-particle states 
should be orthonormal. For this, we need to define the theory in a finite but large region 
whose dimensions are L^, Ly and along the three spatial axes. This gives 

C=^^. (2.38) 

Then 



MX) \n,p,) = \j^j;^f''''Us{y,n,p,) |0) , (2.39) 



where ipu denotes the term in Eq. ()2.36|) that contains the t/-spinors. Similarly, 

[n,p, 



MX) = \^^e^'-'^UMn,n) (0| . (2.40) 



2EnLxLz 
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2.4 Inverse bet a- decay 
2.4.1 Preliminary comments 

Now we have all the tools required for calculating the cross-section of the inverse beta- 
decay process in a magnetic field. We consider the possibility that the neutrons may be 
totally or partially polarized in the magnetic field, and find the cross-section as a function 
of this polarization. The neutrinos are assumed to be strictly standard model neutrinos, 
without any mass and consequent properties. The presence of the magnetic field breaks 
the isotropy of the background, and a careful calculation in this background reveals a 
dependence of the cross-section on the incident neutrino direction with respect to the 
magnetic field. 

Considerable work has been done on the magnetic field dependence of the URCA 
processes which have neutrinos in their final states ^Tj El El ISl HH E] • An angular 
dependence obtained in the differential cross-section of these reactions imply that in a 
star with high magnetic field, neutrinos are created asymmetrically with respect to the 
magnetic field direction. The process that we consider, on the other hand, have neutrinos 
in the initial state. So this process influences the neutrino opacity in a star. 

Some calculations of this process exist in the literature. Roulet [IH]; as well as Lai 
and Qian performed the calculation by assuming that the magnetic field effects enter 
only through the phase space integrals, whereas the matrix element remains unaffected. 
Gvozdev and Ognev considered the final electron to be exclusively in the lowest 
Landau level. Arras and Lai [IH] calculated only the angular asymmetry, and only to the 
first order in the background magnetic field. Some earlier calculations [HTl EHl did not 
take neutron polarization into account. 

In this section the calculation of the cross-section for the inverse beta-decay process 
z/g + ra p + e~ in a background magnetic field is presented in full detail. The calculations 
involve evaluating the matrix element using spinor solutions of the electron in a magnetic 
field, taking all possible final Landau levels into account, including the possibility of 
neutron polarization, and performing the calculations to all orders in the background field 
in the 4-fermi interaction theory. The magnetic field might provide a net polarization of 
the neutrons, which is taken into account. However, the magnitude of the field is assumed 
to be much smaller than rrip/e, so its effects on the proton spinors are ignored. The electron 
spinors, on the other hand, are the ones appropriate for the Landau levels. Thus, we can 
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write the process as 



iy,iq)+niP)-^piP')+eipl,n'). 
2.4.2 The .S-matrix element 



(2.41) 



The charged current interaction Lagrangian for this process is given by Eq. (jl.3|l . 
Using it in first order perturbation, the S-matrix element between the final and the initial 
states of the process in Eq. (|2.4ip is given by 



S 



fi 



V2G 



X (p(P') ii(p)i,{Gv - G^75)V^(n) n{P)) . 



(2.42) 



For the hadronic part, we should use the usual solutions of the Dirac field which are 
normalized within a box of volume V , and this gives 

^i{P'-P)-X 



(2.43) 



using the notations S = Pq and S' = Pq. For the leptonic part, we need to take into 
account the magnetic spinors for the electron. Using Eq. ()2.40j) . we obtain 



(2.44) 



Here Q = qq. The symbol m was used previously to denote the mass of an arbitrary 
charged particle, from now onwards it will denote the mass of an electron in this chapter. 
Putting these back into Eq. ()2.42j) and performing the integrations over all co-ordinates 
except y, we obtain 



Sf, = {2n)%\P + q-P'-p') 



(2.45) 



mV 2£V 2£'V2En'L^L, 
Here, 6^ implies, in accordance with the notation introduced earlier, the 5-function for 
all space-time co-ordinates except y. Contrary to the field-free case, we do not get 4- 
momentum conservation because the y-component of momentum is not a good quantum 
number in this problem. The quantity is the Feynman amplitude, given by 



(2.46) 



= ^2Gp u^pj{P')-f,{Gv - GAl5)u^n){P) / dy e»-^ U(^,){y,n' ,j^)YLu(^.A<l 
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using the shorthand 



'^y ~^ % ^y 



(2.47) 



The transition rate in a large time T is given by \Sfi\^ /T. From Eq. ()2.45j) . using the 
usual rules like 

T 



5{£ + VL-£' - 
S{Px + qx - - p'^ 



2n 



6{S + Q-S' - 



2 



5(p, + g.-p^-p;) =-^5(p, + g,-p;-p;), 

Zn 



we obtain 



, En' + m 



v^nss'En' 



(2.48) 



(2.49) 



2.4.3 The scattering cross-section 

Using unit flux 1/V for the incident particle as usual, we can write the differential 
cross-section as 



da = V 



\S 



fi\ 



T 



-dp, 



where dp, the differential phase space for final particles, is given in our case by 



(2.50) 



(2.51) 



Therefore 



da = V 



T (27r) 



■ dp^ dp. 



= 64^^^^ + ^-^'-^' 
The square of the matrix element is 



(27r)3 

En' + m 



nss'En' 



'^dp'^dp'^d'p'. 



(2.52) 



(2.53) 



where Hny is the hadronic part and l'^^ the leptonic part, whose calculation we outline 



now. 



For the hadronic part, we can use the usual Dirac spinors because of our assumption 
that the magnetic field is much smaller than rrip/e. We will work in the rest frame of 
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the neutron. Due to the presence of the background magnetic field, the neutrons may be 
totally or partially polarized. We define the quantity 

Ar{+) - A^(-) 

S ^ ^- ^ , (2.54) 

where Nj^^ denote the number of neutrons parallel and antiparallel to the magnetic field. 
Then 

H,u = + 5)^i+) + ^(1 - S)H(~J , (2.55) 

where Hj^^ denotes the contribution calculated with spin-up and spin-down neutrons 
respectively. Either of these contributions can be calculated by using the spin projection 
operator, which is |(1±7573) for up and down spins. A straight forward calculation then 
yields 

H,, = 2{Gl + G\){P,Pl + P,P; - g^,P ■ P') 

+2{Gl - G\)mnmpg^, + AiGvGAS^yXpP^P'" 
- S [AGvGAmniP'^g3u + Pu93^. - P-^gf^u) + 2te^,3aR°'] , (2.56) 

where we have introduced the shorthand 

i?" = (G^ + G'^)m„P'° - {Gl - G\)m,P'^ . (2.57) 

We have omitted some terms in the expression for H^j^u that involve spatial components 
of the neutron momentum, with the anticipation that we will perform the calculation in 
the neutron rest frame. 

In the leptonic part we should use the magnetic spinors given in section This 
gives 



j dyj dy, e'-y^y*-y^ lT[Pu{y,y.,n\j^)^^i-i-L\ , (2.58 



where Pu denotes the spinor sum for the electrons, given in Eq. ()2.27p . We now have to 
perform the integrations over y and y^. Each of these variables should be integrated in 
the range —\Ly to +\Ly. However, since we will take the infinite volume limit at the end 



as usual, we let Ly ^ oo and use the result 

r+oo 



/+00 . 
dye'^yin{y)=t''V2^In{x). (2.59) 
-oo 
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This gives 



27r 



where 



eB {En' + m) 



(2.60) 



WeB 
-2V2n'eBg^In' 



pT) + 



^eB, 



In' 



U., 



'eB, 



'eB, 



(2.61) 



Thus, 



fi 



ei3 + m) 



(G^V' + G\){P -AP'-q + P'-AP-q) 

cjn' -r lit-) L 

-(G^ - G^)m„mpg ■ A - 2GvGa{P ■ A P' ■ q - P' ■ A P ■ q) 
+S(2GvGAm^{P' ■ Aq, + P' ■ gA,) -A,q-R + g,A ■ r] 



(2.62) 



We now choose the axes such that the 3-momentum of the incoming neutrino is in the 
x-z plane. We will also assume that |P'| <^ for the range of energies of interest to us. 



In that case, it is easy to see that the terms involving \/2n'eB drop out, and we obtain 

\Gl + 3G^)nAo + {Gl - G\)q,A, 



'"fi 



2 2 27r m„mp 
= 8G^ X 



eB En' + m 



-2GaS({Gv - Ga)^A-, + {Gv + G^)g,Ao)] .(2.63) 



We now put this expression into Eq. (|2.52p and calculate the total cross-section by 
performing the integrations over different final state momenta appearing in that formula. 
First we integrate over P'^ and P'^. These appear only in the momentum conserving 5- 
function. Integration over them therefore just gets rid of the corresponding 5-functions. 
For the integration over p^., we refer to Eq. (j2.13p . Since the center of the oscillator 
has to lie between —\Ly and \Ly, we conclude that —\LyeB < p'^ < \LyeB. Thus the 
integration over p'^ gives a factor LyeB. 

Putting back into Eq. (j2.52p and using V = L^LyLz-, we obtain 

Gl5{D + n- En') 



da 



An 



VtEn' 



{Gl + ?>Gl)nAo + {Gl - G\)q,A, 

+2GaS({Gv - Ga)^A, + {Gv + GA)g.Ao)] dP'ydp', , 

(2.64) 
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where D is the neutron-proton mass difference, m„ — m^. 

We next perform the integration over Py. In the integrand, it occurs only as the 
argument of the functions and In-i- The functions are orthogonal in the sense that 

n+OO 



(2.65) 



/ + 00 
da In{o)Ini{a) = y eB 5n 



This property can be used to perform the integration over Py. We have already remarked 
that the term proportional to \/2n'eB in Eq. (|2.6ip does not contribute. From other two 
terms, we obtain 



eB 
eB 



9n'P'\\ + K'flP\ 



where 



2 — 5r,l 



n',0 



(2.66) 



(2.67) 



gives the degeneracy of the Landau level. Notice the appearance of the Kronecker delta, 
5n'fli in the expression of Eq. ()2.66|) . The reason for this is that, while two terms of Eq. 
(j2.61|) contribute in the integral for n' 7^ 0, only one of them contributes for n' = since 
/_i = 0. 

The final integration is over p'^. Writing the argument of the remaining (5-function in 
terms of p'^, we find that the zeros occur when 



p'^ = p'^ = ±J(D + ny - m2 - 2n'eB . 



(2.68) 



D + Q 



= {6{p',-p'^) + 6{p',-p'_)). (2.69) 



Therefore, 

6{D + n- En') 

^{D + - m2 - 2n'eB 

In the integration, the terms proportional to p'^ in the integrand receive equal and opposite 
contributions from the two 5 functions and cancel. For the other terms, independent of 
both the contributions are equal. So we obtain 

eBGl 



271 



gn'[{Gl + ?,G\) + 2GaS{Gv + Ga) cosO] 
+5n',o[{Gl - G\) cose + 2GaS{Gv - Ga)} 



D + Q 



{D + fi)2 - m2 - 2n'eB 

(2.70) 
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where we have defined the direction of the incoming neutrino by the angle 6, with 



q, = ncos9. (2.71) 

In Eq. ()2.7()j) . we have denoted the cross-section by an' because the electron ends up in 
a specific Landau level n'. The total cross-section is then given as a sum over all possible 
values of n', i.e., 

^ = E ^n' = E 9n'{{Gl + 3Gl) + 2GaS{Gv + Ga) cosO} 

n'=0 n'=0 

+Sn',o{iG'v - Gl) cose + 2GaS{Gv - Ga)}] , • 
^ ^{D + Q)^ - m2 - 2n'eB 

(2.72) 

The possible allowed Landau level has a maximum, TT-'^ax' which is given by the fact 
that the quantity under the square root sign in the denominator of Eq. ()2.72|) must be 
non- negative, i.e., 

/ 



= int 



max 



^ {D + nf-m^\. (2.73) 



\2eB 

Eq. ()2.72p gives our result for the cross-section of the inverse beta decay process. Some 
properties of this formula are worth noting. 

For unpolarized neutrons, S = 0, the cross-section for n' 7^ does not depend on the 
direction of the incoming neutrino. The same is not true if the electron ends up in the 
lowest Landau level. The cross- section will be asymmetric in this case. 

All terms in the cross-section which depend on 5* have a common factor Ga- The 
reason is that, if Ga were equal to zero, the interaction in the hadronic sector would have 
been spin-independent. 

If the final electron is in the lowest Landau level and the initial neutrino momentum 
is antiparallel to the magnetic field, Eq. ()2.72p shows that 



2tt 



(2.74) 



Note that the vector coupling Gy does not contribute to the cross- section in this limit. 
This can be understood easily. The neutrino spin is along the +z direction whereas the 
electron spin in the lowest Landau level must be in the —z direction. Thus there is a 
spin-fiip in the leptonic sector. Conservation of angular momentum then implies that 
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Figure 2.1: Enhancement of cross-section in a magnetic field for an initial neutrino energy of 
lOMeV. Different panels show the results for different net polarizations of the neutrons. The solid 
and the dashed lines correspond to the initial neutrino momentum antiparallel and parallel to the 
magnetic field. Each curve has been normalized by the cross-section in field-free case for the same 
values of S and cos6'. The horizontal dashed lines represent unity in the vertical scale. 



there must be a spin-flip in the hadronic sector as well. In the non-relativistic limit for 
hadrons that we have employed, this can occur only through the axial coupling. 

If further we consider totally polarized neutrons, i.e., S* = 1, we see that o"o vanishes. 
Again, this is a direct consequence of angular momentum conservation. Since both initial 
particles have spin up, angular momentum conservation requires both final particles in 
spin up states as well. But the spin-up state is not available for the electron in the lowest 
Landau level. 

It is instructive to check that the result obtained in Eq. (|2.72|) reduces to the known 
result for the field-free case. The contribution specific to the zeroth Landau level vanishes 
in the limit B owing to the overall factor of eB. The other terms also have the factor 
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eB, but in this case we also need to sum over infinitely many states. This gives 
a = UgI. + 3Gl) + 2GaS{Gv + Ga) cos 6 

d + vl d + vl 




{D + fi)2 - m2 - 2n'eB 2J{D + - m 



,2 



(2.75) 



For i3 — s> 0, the last term vanishes, and we can identify n'^^^ as the integer for which the 
denominator of the summand vanishes. Thus we obtain 

eBGl r , „o „o . „ „ . „i /""max . , D + il 



TT 



r<2 



{GI + 3G\) + 2GaS{Gv + G a) cose] I 



dm! 



{D + ny -m? - 2n'eB 



{G'y + SGl) + 2GaS{Gv + Ga) cosO {D + n)J{D + n)^ - , (2.76) 

7]- L J V 

which is the correct result in the field-free case. 

In Fig. 12.11 we have plotted the ratio of the cross-section to its corresponding value 
at i3 = as a function of the magnetic field. The plots have been done for unpolarized 
{S = 0) as well as totally polarized neutrons along {S = 1) and opposite {S = —1) to 
the magnetic field, with the initial neutrino momentum parallel and antiparallel to the 
magnetic field. For 5* = — 1, we find that neutrinos parallel to the magnetic field have a 
smaller cross-section than those antiparallel to the field, and the difference is pronounced 
for large fields. For S = 0, the situation is just reversed. For 5* = 1, if the magnetic 
field is high enough so that ^'max = 0, we see that the cross-section vanishes for neutrino 
momentum antiparallel to the field. The reason for this has already been discussed. 

For the inverse beta decay process, which is symbolically written as in Eq. ()2.4H) . the 
cross-section has been calculated earlier by several authors [IHl EZ] ■ They assumed that 
the matrix element remains unaffected by the magnetic field, only the modified phase 
space integral makes the difference in the cross- sect ion. The magnetic field effects enters 
only through the following modification of the phase space integral and the spin factor of 
the electron: 

d^p eB 



T.9n'Jdp,, (2.77) 



(27r)3 (27r)2 ^ 

where Qn' is the degeneracy of the n'-th Landau level as explained in Eq. ()2.67|) . The sum 
over n' is restricted to the region 

, (D + Q)'^~mj 
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where Q, is the neutrino energy. The results they obtained is the same as the term 
proportional to Gy + 3G\ that we obtained. 

Subsequent calculations incorporated the modification of wave functions. Arras and 
Lai [45 , while still treating the nucleons as non-relativistic, used the non-relativistic 
Landau levels as well as the finiteness of the recoil energy for the proton. They found the 
cross-section and went on to derive expressions for the neutrino opacity. From the final 
expressions, one can only recognize the terms linear in B. The opacity was calculated also 
by Chandra, Goyal and Goswami Like the previous authors, they also considered the 
contribution to the opacity from other reactions like neutrino-nucleon elastic scattering. 

2.4.4 Consequences of neutrino energy spread 

The enhancement factor in Fig. 12.11 shows some spikes. They appear at values of the 
magnetic field for which the denominator of Eq. (j2.72j) vanishes for some n'. For field 
values larger than this, that particular Landau level does not contribute to the cross- 
section. To the right of the final spike that appears in the figure, only the zeroth Landau 
level contributes. In other words, the final electron can go only to the lowest Landau level 
for such high values of the magnetic field. The exact value of B for which this occurs 
depends of course on the energy of the initial neutrino. 

We need to make an important point about these spikes. Each spike in fact go all the 
way up to infinity. The finite height of a spike in the figure is an artifact of the finite step 
size taken in plotting it. 

In reality, of course, a cross-section cannot be infinite. In the present case, the cross- 
section remains finite due to the fact that the initial particles cannot be exactly monochro- 
matic due to the uncertainty relation. For an example if we concentrate on the initial 
neutrino alone then there must be a spread in its energy, which can be represented by a 
probability distribution defined by 



In that case, the cross-section in a real experiment should be written in the form 




(2.79) 




(2.80) 



where a{Q) is the expression derived in Eq. ()2.72|) for a single value of energy. 
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Figure 2.2: Cross-section for unpolarized neutrons as a function of the magnetic field for a flat energy 
distribution, normalized to the cross-section in the field-free case. The initial neutrino momentum is 
along the magnetic field, and energy is 10 MeV. The energy spread — is 0.2 MeV for the left 
panel and 2 MeV for the right panel. 



As an illustration, we consider the case of unpolarized neutrons (5 = 0), and take a 
flat probability distribution of initial neutrino energy, viz.. 



if fii < < 1^2, 



ih-ih 

otherwise. 
Then the integration of Eq. ()2.8()j] gives 



(2.81) 



a 



eBG} 



27^(^]2 - ^^1 



F(^]2)-m) 



(2.82) 



where 

n' 

' "max 

m) = E 



n'=Q 



X J{D + ny -m"^ -2n'el3, 



(2.83) 



with n'^^^ determined by Eq. ()2.73|) . Fig. 12.21 shows the variation of this quantity with 
the magnetic field for cos^ = 1. In this figure, we normalize the cross-section by ;B = 
cross-section with the energy distribution of Eq. (j2.81|) . which is 



a{B = 0) 



G'b{Gv + SG^) 
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(2.84) 



3TT{n2 - ill) 

Keeping the central value of neutrino energy as 10 MeV as before, we have drawn these 
plots for two different values of the spread, as mentioned in the caption. For the smaller 
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value of the spread in particular, the graph looks very similar to that drawn in Fig. 12. H 
but the difference is that now the height of the spikes denote the actual enhancement, and 
is not an artifact of the plotting procedure. For the higher value of the energy spread, we 
see that the spikes have smoothened out. 

2.5 Asymmetric emission of neutrinos from a proto- neutron star 

From what has been discussed so far we know that the inverse beta decay cross-section 
is direction sensitive through its dependence on the angle 9, the angle between the neutrino 
3-momentum and the magnetic field direction. This anisotropic effect in the cross-section 
can have interesting astrophysical applications in explaining the high velocities of pulsars, 
of the order of 450 ±90 Kms~^. Typical pulsars have masses between I.OMq and 1.5Mq, 
i.e., about 2 x lO^^g. The momentum associated with the proper motion of a pulsar would 
therefore be of order 10^^ gcm/s. On the other hand the energy carried off by neutrinos 
in a supernova explosion is about 3 x lO^^erg, which corresponds to a sum of magnitudes 
of neutrino momenta of lO^^gcm/s. Thus an asjmametry of order of 1% in the distribution 
of the outgoing neutrinos would explain the kick of the pulsars. It has been argued that 
an asymmetry of this order in the distribution of outgoing neutrinos can be generated by 
the anisotropic cross-sections of the various neutrino related processes in presence of a 
constant magnetic field [ini EH IS21 UHl HSl IHIj • A simple calculation can make the point 
clear. 

The typical size of a neutron star is about R 10 km. Thus, when a neutrino, 
generated at the core reaches a density where its mean free path is about -R, it escapes 
from the star. Therefore the condition for the neutrino to escape can be written as 



where is the neutron number density. We already observed that a is direction depen- 
dent. Therefore, the value of n„ on the "neutrino sphere" depends on the direction as 
well, and the surface is no longer a sphere. Different values of will correspond to dif- 
ferent temperatures. Thus, neutrinos will be emitted with different momenta in different 
directions. This can result in a kick to the star. 

To estimate the magnitude of the kick, let us abbreviate Eq. ()2.72|) as 



1 



(2.85) 




a = eBGl{a + bcos9) . 



(2.86) 
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If we now consider the directions ^ = and 6 = tc, the difference in neutron density on 
the corresponding points on the neutrino surface is given by 

_ 1 / 1 1 A _ 26 

neglecting corrections of order b/a. 

The neutron gas in a typical proto-neutron star can be considered to be non-relativistic 
and degenerate. The number density of neutrons is thus given by jHS] 



1 H h 



(2.88) 



where is the Fermi momentum, and we have neglected higher order terms in the 
temperature. This gives 



dn„ mi ( T 



n 



(2.89) 



dT 3 \pf ) 

So the temperature difference between the points on the neutrino surface in the ^ = 



and 6 = 7T directions is 



^^ = 4rs|k (2,90) 



The momentum asymmetry can now be written as 

A\q\ 1 4AT 



(2.91) 



\q\ 6 T ' 

where we have assumed a black body radiation luminosity (oc T^) for the effective neutrino 
surface. The factor 1/6 comes in because the asymmetry pertains only to z/g, whereas 6 
types of neutrinos and antineutrinos contribute to the energy emitted. This gives 

Alqrl A pf b , , 
Lli = £_L 2 921 

|g| mlT^eBGlo?R' ^ ' 

To find pf, we use the leading term in Eq. ()2.88p and estimate n„ from the equation 

P(l -Ye) , , 

rin = , 2.93 

rrip 

where Ye is the electron fraction and p is the mass density. Taking = 1/10, we obtain 

Pf = 24p}(^ MeV . (2.94) 
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where pu is the mass density in units of lO^^gcm ^. Putting this back in Eq. ()2.92|) . we 
obtain 



27nV3 g-i'7.-2 \b[ 



(2.95) 



where Bu = 0/(10^'' Gauss) and Tmcv = T/{1 MeV). For i7 ^ me which is the relevant 
case and S* = as a typical value of neutron polarization we can find out the values of 
the constants a and b. The values of both of them will in general depend upon n^^x but 
if we take only the contribution from the n' = level then 
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using Gy = 1 and Ga 



2-K 

-1.26. This gives 
A\q 
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(2.96) 



(2.97) 



(2.9J 



Obviously, with reasonable choices of p, B and T, it is possible to obtain a fractional 
momentum imbalance of the order of 1% which is necessary for explaining the pulsar 
kicks. 



2.6 Conclusion 

This chapter started with a discussion on neutrino-electron scattering and other pro- 
cesses related to it by crossing. The basic theory of spinor solutions in an external mag- 
netic field has been exposed with some detail in next two sections. A formal field theo- 
retical framework has been presented based on which cross-sections of any process, which 
contains charged particles in the initial or the final state, can be calculated. The main 
calculations concerning the inverse beta decay cross-section start from the fourth section. 

The calculations show that for the inverse beta-decay process, even for unpolarized 
neutrons, the cross-section depends on the direction of the neutrino momentum. This 
asymmetry is not surprising since the background magnetic field makes it an anisotropic 
problem. A similar asymmetry has been noted for URCA processes PE], where the neu- 
trino is in the final state. From Eq. ()2.72|) it is seen that the anisotropy in the cross-section 
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comes only from the zeroth Landau level contribution of the electron 03]. However, for 
the higher levels, there is a cancellation between the two possible states in a Landau level 
which washes out all angular dependence in these levels, provided the neutrons are un- 
polarized. The asymmetry in cross-section will therefore come only from the the zeroth 
Landau level state of the electrons and its amount will depend on the relative contribu- 
tion of this state to the total cross-section. If the magnetic field is so high that only the 
zeroth Landau level state can be obtained for the electron, the asymmetry will be large, 
about 18%. For smaller and smaller magnetic fields, the asymmetry decreases with new 
Landau levels contributing. For polarized neutrons, however, there is an asymmetry even 
in the field-free case as is evident from Eq. (I2.76|) . In presence of a magnetic field, the 
asymmetry will in general depend on the magnetic field, as it appears from various plots 
in Fig. O 

From Eq. (j2.76|) it is seen that the inverse beta decay cross-section has a smooth B 
limit. From the curves in Fig. 12. H it is observed that the cross-section has spikes, which go 
all the way to infinity, for specific values of the magnetic field magnitude. Subsection l2.4.4l 
deals with the issue of these infinite spikes. It is shown that if instead of a monochromatic 
neutrino the initial neutrino has an energy spread, as is expected in a real situation, then 
the spikes gets smoothed up. This fact is evident from Fig. 12.21 

The anisotropic nature of the inverse beta decay cross-section can have far reaching 
consequences for neutrino emission from a proto-neutron star as has been shown in the 
section concerning asymmetric emission of neutrinos. It has been discussed in the lit- 
erature that the presence of asymmetric magnetic fields can cause asymmetric neutrino 
emission from a proto-neutron star [53]. However, the calculations in this chapter show 
that even with a uniform magnetic field, neutrino emission would be asymmetric because 
of the cos ^-dependent terms in the cross-section. 
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CHAPTER 3 



Observable effects due to 
intermediate virtual charged 
particles in a magnetic field 

3.1 Introduction 

In quantum field theoretic calculations, the spinors given in the previous chapter in 
Eq. ()2.24j) and Eq. ()2.2(jj) should be used if the charged particle appears in the initial or the 
final state of a physical process. If, on the other hand, the charged particle appears in the 
internal lines of the Feynman diagram of the process then we should use its propagator. 

This chapter is on propagators for charged particles, mainly fermions, in presence of 
a background uniform magnetic field. The calculations will be extrapolated to the case 
where there is a background medium also. For the medium modification of the propagators 
a brief discussion on statistical field theory (commonly called finite temperature field 
theory) is discussed in the appendix O The Schwinger propagator is introduced in this 
chapter and is followed by a discussion on the phase factor accompanying it. A detailed 
discussion on the derivation of the Schwinger propagator is supplied in The phase 
factor of the Schwinger propagator is dealt in some detail. In most of the works this phase 
factor is treated with irrelevance but the reason is rarely explained. 

In the chapter^ we noticed that the standard model neutrinos do not have any mag- 
netic moment and so they cannot interact with external magnetic fields. In going beyond 
standard model i.e. by including right handed neutrinos we can generate a magnetic mo- 
ment of the neutrinos. Magnetic fields can also affect the neutrino properties in a different 
way i.e., through the charged particle propagators as discussed in this chapter. This in- 
direct effect of the magnetic fields from the propagators can affect neutrino self-energy 
calculations. Without doing actual loop calculations it is possible to predict the form of 
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the self-energy expressions in a magnetic field and in a magnetized medium using a purely 
general form-factor analysis. Although such an analysis is purely formal but still it can 
predict interesting effects, like the dependence of the neutrino self-energy on the angle 
between the neutrino propagation direction and the external magnetic field direction, in 
a magnetized medium. The fact that the neutrino self-energy is sensitive to the angle 
between the neutrino propagation direction and the magnetic field direction can have far 
reaching consequences. In the field of neutrino oscillations the angular dependence of the 
self-energy can modify the resonant level crossing condition of the neutrinos. The res- 
onant condition itself becomes sensitive to the angle between the propagating neutrinos 
and the external magnetic field which can have an interesting astrophysical application 
like the asymmetric emission of neutrinos from a newly formed neutron star. 

In the two following sections the fermion propagator in a magnetic field in vacuum and 
in a thermal medium will be discussed. In the remaining sections the topic of neutrino 
self-energy in a magnetic field and the effect of magnetic fields on neutrino oscillations will 
be discussed briefly. The chapter ends with a conclusion which summarizes the various 
things discussed in this chapter. 

3.2 Propagators in vacuum 
3.2.1 Furry propagator 

There are two ways to write the propagator. The first is to start with the fermion field 
operator ipi^) written in terms of the spinor solutions and the creation and annihilation 
operators, and construct the time ordered product, as is usually done for finding the 
propagator of a free fermion field in the vacuum. The algebra is straight forward and 
yields the result 



where is the positive root obtained from Eq. ()2.16|) . The notation requires some 
explanation. Instead of x use has been made of X as the coordinates in the left hand 
side of the above equation. This style of writing is carried over from section 12.21 It is 
motivated by the fact that in calculations where we employ these propagators we often 
have to write down p ■ = p^x + pzZ in the explicit form where the quantity x appears 






s 
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as a coordinate and not as the 4-vector itself. The notation S^{X,X') specifies that the 
propagator is in vacuum and in presence of a uniform background magnetic field whose 
magnitude is B. The spin sum appearing in Eq. (j3.ip is given in Eq. (|2.27|) . The resulting 
propagator is called the propagator in the Furry picture. 

3.2.2 Schwinger propagator 

Alternatively, one uses the propagators introduced by Schwinger [SB]- Schwinger's 
calculation of the fermion propagator relies on a functional procedure and it is written in 
the form 

tSUx. x') = nix, x') J ^e-^-(^-^')^5]^(p) , (3.2) 

here x stands for the coordinate 4-vector as usual. S]^{p) is expressed as an integral over 
a variable s, usually (though confusingly) called the 'proper time': 

iSl{p)= dse^^^'^Gip.s). (3.3) 
Jo 

The quantities $(p, s) and G{p, s) can be written in the following way : 



<l>(p, s) = is 



, tan(ei3.)^^_^A_^l^l^ 



eBs 



= {1 + ia,tan{eBs)){^ii + m) + sec'^{eBs)^± . (3.5) 

In the above expressions az = ^7^7^. In a typical loop diagram, one therefore will have 
to perform not only integrations over the loop momenta, but also over the proper time 
variables. iSsip) is manifestly translation invariant. Most of the calculations involving 
the Schwinger propagator use this translation invariant part only, while the other term 
k{x^x') appearing in the propagator in Eq. ()3.2|) remains irrelevant. 

Not going into any detailed description of k(x, x') we can simply understand its neces- 
sity in the propagator. A propagator connects two points in space-time. In presence of 
a background gauge field the gauge transformation property of the fields of the charged 
fermions are different at two different space-time points. Unless there is some factor in the 
propagator which can connect these two fields at different space-time points with different 
gauge transformation properties, the calculations involving charged fermion propagators 
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will not be manifestly gauge invariant. In the limit the background magnetic field goes 
to zero Eq. ()3.3|) shows that iS^{p) reduces to the vacuum fermionic propagator. This 
indicates that iS]^{p) in Eq. (|3.2|) do not carry any information about gauge invariance 
of the background gauge field, it is only the translation invariant modified version of 
the vacuum fermionic propagator. Therefore the gauge transformation property of the 
propagator must be related to k{x,x'), and its presence is necessary for the overall gauge 
invariance of physical processes. The gauge transformed fields comes with phase factors 
where the phase depends upon the space-time point where the gauge transformation is 
made. The fermionic fields at two different space-time points will therefore have two 
different phase factors. To make a connection between them k,{x, x') must also be some 
form of a phase. Conventionally it is named the phase-factor. 
The phase factor is generally written as 

k{x, x') = exp {iel{x, x')} (3.6) 

where 



j(x, x') = r de k(o + If^a^ - x' 



(3.7) 



where is the background gauge field, and F^,^ is the field strength tensor. From Eq. ()3.7|) 
we notice that the phase factor breaks the translation invariance of the propagator. 
For a constant background field we can always write the gauge field as 

A,iO = -lF,,C + d,m, (3.8) 

where A(^) is an arbitrary well behaved function and depends upon our choice of gauge. 
Using the above relation in conjunction with Eq. (|3.7|) we can simplify the integration 
appearing in the phase factor as 



i{x,x')= rde I-If.^x"' + d.xio 

Jx' L ^ 



(3.9) 



Using the constancy of the field strength tensor the above expression can be written as 

I{x, x') = ^x'T^^x" + X{x) - A(x') . (3.10) 

From Eq. ()3.10|) we can immediately see if we set x = x', in other words if we integrate 
over any closed contour in space-time I{x,x') vanishes. Thus I{x,x') connecting two 
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points in space-time is independent of the path joining them, and as a result the phase 
factor of the Schwinger propagator joining the points x' and x in Eq. ()3.6p is also path 
independent. 

Utilizing the path independence of the phase factor of the propagator the general 
convention is to choose a straight line path connecting the two points x' and x. Points on 
this path are represented by 

e = (l-cy'^ + ca:^ (3.11) 

where the parameter ranges from to 1. Using Eq. (j3.7p and the straight line path 
given in Eq. (j3.11|) we get 

J(x,x') = rd^A^iO + i f' dCix^ - x'^F^^ix" - x"") , (3.12) 

Jx' I Jo 

Jx' 

Now using Eq. ()3.8|1 for the gauge field we can retrieve Eq. ()3.1()j) . 

From Eq. ()3.1()|1 it is clear that the phase factor is dependent on the form of the 
function A(^), that is to say the fermion propagator is dependent on the gauge in which 
the constant background magnetic field is specified. Suppose we are working in such a 
gauge that A(^) = 0, and then we make a gauge transformation of the background field 
as 

A^^A^ + d^m (3.13) 
then from Eq. (|3.2j) it follows that the fermion propagator will transform as 

iS^{x, x') — > exp{ieX{x))iS^{x, x') exp{—ieX{x')) . (3-14) 

under the gauge transformation. As because we are working in presence of a background 
gauge field, the fields of the charged particles and their propagators both become back- 
ground gauge dependent. This is the reason why the phase factor arises in the expression 
of the propagator. A detailed study of the phase factor appears in appendix IB. 11 

3.3 Schwinger propagator in presence of a medium 

In presence of a medium the propagators of the elementary particles get changed. 
There are various formulations of statistical field theory, like the imaginary time method, 
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the real time formalism, to name the most prominent ones. In the present work we will 
be using the real time formalism developed in the canonical approach [22]. Some of the 
important points of the above formalism is discussed in appendix O 

In this formalism it is not difficult to obtain the effects of a magnetized medium on the 
Schwinger propagator. In the real-time formalism, the 1-1 component of the propagator 
iS{p) involving the time-ordered product'^ can be written in terms of the free propagator 
iSo{p): 

iS{p) = iSoip) - r]Fip ■ u) [iSoip) - iSo{p)] , (3.15) 

where 

So{p) =loSl{pho, (3.16) 

and rip{p ■ u) contains the distribution function for particles and antiparticles: 

Vrip ■ u) = Q{p ■ u)fFip, /X, (3) + 0(-p ■ u)fF{-p, -/i, P) ■ (3.17) 

Here, G is the step function which takes the value +1 for positive values of its argument 
and vanishes for negative values of the argument, is the 4- vector denoting the center- 
of-mass velocity of the background plasma. Conventionally u is normalized in such a way 
that in the rest frame of the medium 

u^ = {l,0). (3.18) 
fp denotes the Fermi-Dirac distribution function: 

/Hp,/x,/3) = ^,(^J^)^^ . (3.19) 
In a similar manner, the 1-1 component of the propagator in a magnetized plasma is given 

by jsaisHi 

iSsip) = iS^ip) - r]F{p ■ u) iSsip) - iSsip) ■ (3.20) 

In the Schwinger proper-time representation, this can also be written as an integral over 
the proper-time variable s: 

iSB{p)= dse'^^P^'^G{p,s)-7]F{p-u) ds e'^^P''^G{p,s) , (3.21) 

Jo J-oo 



^It should be mentioned here that other orderings also appear in the evaluation of general Green's 
functions. The other orderings give rise to other propagators as discussed in appendix lUl We will not 
need those other propagators here. 
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(a) (b) 

Figure 3.1: One-loop diagrams for neutrino self-energy in a magnetized medium. Diagram 
b is absent if the background contains only a magnetic field but no matter. For legends and 
related diagrams, see the caption of Fig. 11.11 



where s) and G{p,s) are given by the expressions in Eq. ()3.4p and Eq. ()3.5p . Sym- 
bohcally the above equation can be written as 

^SB{p)=^Sl{p) + Sl{p) (3.22) 

where iS^{p) is the vacuum propagator in a uniform magnetic field and is given by 
Eq. (|3.3p . The other part S^{p) of the propagator carries the medium effect, and is given 
by 

/oo 
ds e*(P'")G'(p, s) . (3.23) 
-oo 

It is straight forward to see that when B = 0, the propagator in Eq. ()3.3p reduces to 

poo p . , -, 

'i'Soip) = J ds exp is [p^ — m"^ + iej (^ -|- m) 



p2 _ ^2 _j_ 



(3.24) 



which is the normal vacuum Feynman propagator. In the same limit, the background 
dependent part reduces to 

S^{p) = -2^ 5{p^ - m^)T]F{p ■ u){i) + m) , (3.25) 

which is what we expect in a medium in absence of a magnetic field. 

3.4 Neutrino self-energy 



In chapter ^ the self-energy of the neutrino was discussed in a background magnetic 
field and an expression of the self-energy was supplied in Eq. ()1.2|) . Here we take up 
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the issue and briefly discuss about the form of the self-energy expression [21]. Skipping 
detailed calculations using the Schwinger propagators a form-factor analysis of the neu- 
trino self-energy is given in this section and its implications are highlighted in the case of 
neutrino oscillations. 

It is already mentioned in chapter ^ that the simplest physical quantity where back- 
ground magnetic field effects appear through virtual lines of charged particles is the self- 
energy of the neutrino. The 1-loop diagram for the self-energy is given in Fig. 13.11 In the 
present discussion we take the case of neutrinos in a medium with a magnetic field. From 
this general approach we will be able to find the form of the self-energy of the neutrino 
in absence of a medium also. 

It is easy to see how the self-energy might be modified within a magnetized plasma. 
In the vacuum, the self-energy of a fermion has the general structure 

= a7^g^ + 6, (3.26) 

which is the most general form dictated by Lorentz covariance. Here, a and h are Lorentz 
invariant, and can therefore depend only on q^. In the presence of a homogeneous medium, 
the self-energy will involve the 4- vector u^. Further, if the medium contains a background 
magnetic field, the background field F^j_i, also enters the general expression for the self- 
energy. These new objects, and F^^, enter in two different ways. First, any form factor 
now can depend on more Lorentz invariants which are present in the problem. Second, 
the number of form factors also increases, since it is possible to write some more Lorentz 
covariant terms using and F^y. There will in fact be a lot of form factors in the most 
general case. However, if we have chiral neutrinos as in the standard electroweak theory, 
the expression is not very complicated 



^B{q) = {aiq^ + Wu^ + a^q'^F^, + b^u'' F^, + a^q'^F^, + b^u" F^,)-f^'L , (3.27) 

where L is the left-chiral projection operator, and 

F^.u = le^uXpF^r (3.28) 

We first consider the self-energy when the background consists of a pure magnetic 
field, without any matter. Then all 6-type form-factors disappear from the self-energy. 
The dispersion relation of neutrinos can then be obtained by the zeros of ^ — S^, which 
gives 

;i - ai)q>' - aag^F^. - a^q'F^,] ' = . (3.29) 
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It is interesting to note that the terms hnear in the background field all vanish due to 
the antisymmetry of the field tensor. Moreover, the 0203 term is also zero for a purely 
magnetic field. This gives 

(1 - aifq^q^ + alq'qxF^.F''^ + a^qxP^^F''^ + 2a2a3q''qxF^,F^''^ = . (3.30) 

The remaining terms can be most easily understood if we take the 2;-axis along the 
direction of the magnetic field. Then the only non-zero components of the tensor F^^, and 
Ff^i, are given by 

F12 = -F21 = B, Fo3 = -F30 = B , (3.31) 
where we have adopted the convention 

£0123 = +1 • (3.32) 

Thus 

q'^qxF.^F'^' = -{ql + q^B^ = -qlB' , 

q^qxF,^^' = - ql)B' = qp' , (3.33) 

where the notations for parallel and perpendicular products were introduced in the chap- 
ter n The form factor ai can be set equal to zero by a choice of the renormalization 
prescription. So the dispersion relation is now a solution of the equation 

q' - alqlB' + alqp' = , (3.34) 

which can also be written as 

Of course, this should not be taken as the solution for the neutrino energy, because the 
right hand side contains form factors which, in general, are functions of the energy and 
other things. But at least it shows that in the limit B ^ 0, the vacuum dispersion relation 
is recovered. If we retain the lowest order corrections in B, we can treat the form-factors 
to be independent of B and can write [3^ 

n' = q' + (al - al)B\l . (3.36) 
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Calculation of this self-energy was performed by Erdas and Feldman jTU] using the 
Schwinger propagator, where they also incorporated the modification of the TV-propagator 
due to the magnetic field. Importantly, the VT-propagator contains the same phase factor 
as discussed in appendix IB. II Therefore, the phase factors from the charged lepton and 
the IV-lines are of the form k,{x, x')k,{x', x). From appendix IB. II it is easy to see that this 
is equal to unity, and therefore the phase factors do not contribute to the final expression. 
Detailed calculations show that pii] 



where m is the mass of the charged lepton in the internal line. The above equation in 
conjunction with Eq. ()3.36|) gives us Eq. ()1.2j) appearing in the chapter ^ 

Let us next concentrate on the terms which can occur only in a magnetized medium. 
In other words, we select out the terms which cannot occur if the neutrino propagates in a 
background of pure magnetic field without any material medium. This means that, apart 
from the term ^ which occurs also in the vacuum, we look for the terms which contain 
both u'^ and F^^,. Further, if the background field is purely magnetic in the rest frame of 
the medium, u'^F^,^ = since u has only the time component whereas the only non-zero 
components of Fj^^, are spatial. Thus we are left with ISHl '■ 



Once again, setting oi = through a renormalization prescription, we can find the dis- 
persion relation of the neutrinos in the form [3^ : 



where q is the unit vector along q, and we have kept only the linear correction in the 
magnetic field. This form for the dispersion relation was first arrived at by D'Olivo, 
Nieves and Pal (DNP) 59J who essentially performed a calculation to the first order in 
the external field. As for the form factors, bi was known previously, obtained from the 
analysis of neutrino propagation in isotropic matter, i.e., without any magnetic field. The 
result was 133 EOl EH 




(3.37) 




(3.38) 



n= q-hB +bi^\q 



(3.39) 



bi = y/2GF{ne - He) X {ye + pCy) , 



(3.40) 



42 



where 



P = ,.2 2n ^ (3-41) 

M| cos^ Ow 

He, n-e are the densities of electrons and positrons in the medium, 

^^ = |o for.^.„ ^3-^2) 

and cy is defined through the couphng of the electron to the Z-boson, whose Feynman 
rule is 

Z COS 

In other words, in the standard model 

cv = -^ + 2smHw , ^^ = ~l- (3.44) 

The contribution to 63 from background electrons and positrons was calculated by 
DNP ISni- They obtained^ 

63 = -2^2eGp I ^(/^ _ X [y^ + PC J,) , (3.45) 

where /e and /e are the Fermi distribution functions for electrons and positrons, and 



E= Jp^ + ml. (3.46) 



Later authors have improved on this result in two different ways. Some authors |^ have 
included the contributions coming from nucleons in the background. Some others [HZl IS^ 
have used the Schwinger propagator and extended the results to all orders in the magnetic 
field. 

3.5 The role of neutrino oscillations in a magnetized medium 

Neutrino oscillation phenomenology in a medium has produced the important MSW 
effect which is now widely applied in the solar neutrino problem and also in the 

case of neutrinos coming from supernovas. But in most of the cases there is a background 



■*The authors of [S^l used a convention in which e < 0. Here we present the result in the convention 
e > 0. 
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magnetic field and so the following discussion is an important part of neutrino oscillation 
phenomenology. 

Calculation of neutrino self-energy has a direct consequence on neutrino mixing and 
oscillations. Of course neutrino oscillations require neutrino mixing and therefore neutrino 
mass. For the sake of simplicity, we discuss mixing between two neutrinos which we will 
call Vf. and v^. The eigenstates will in general be called vx and v^^ which are given by 



v\ \ _ ( cos9 - sin6' \ 
z/9/ I sin^^ cos6' \iy, 



(3.47) 



We will denote the masses of the eigenstates by mi and 1712, and assume that the neutrinos 
are ultra-relativistic. Then in the vacuum, the evolution equation for a beam of neutrinos 
will be given by 

4M=-M^hV (3.48) 

where the matrix is given by 

2 _ /-|Am2cos2^^ iAm2sin2^\ 
~l^iAm2sin2^ iAm2cos2^j' ^"^^^^ 

where Am^ = m\ — m\. In writing this matrix, we have ignored all terms which are 
multiples of the unit matrix, which affect the propagation only by a phase which is 
common for all the states. 

In a non-trivial background, the dispersion relations of the neutrinos change, as dis- 
cussed in Sec. 13.41 This adds new terms to the diagonal elements of the effective Hamil- 
tonian in the flavor basis, which we denote by the symbol A. As a result, the matrix 
should now be replaced by 

~2 _ / Am2 cos 26 + A,^ i Am^ sin 26 \ 

-\ lAm2sin2^ ^Am^ cos2^ + ' ^^"^^^ 

where the extra contributions are in general different for z/g and v^. The eigenstates and 
eigenvalues change because of these new contributions. For example, the mixing angle 
now becomes 6', given by 

~ Am^sin26' 
tan2e = — — — . 3.51 

Am2 cos 26 + A^^ - A^^ ^ ' 

In a pure magnetic field, the self-energies were shown in Eq. ()3.36|) and Eq. ()3.37|) . The 
quantity m appearing in Eq. ()3.37|) is the mass of the charged lepton in the loop. Thus, for 
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Vei it is the electron mass whereas for z/^, it is the muon mass. Thus A^,^ ^ Ay^. However, 
the difference appears in logarithmic form, and is presumably not very significant. 

In a magnetized medium, however, the situation changes. The reason is that the 
medium contains electrons but not muons. Accordingly, the quantities Ay^ and Ay^ can 
be very different, as seen by the presence of the term ?/e in E<1- (!3.45p . If we take self-energy 
corrections only up to linear order in as done in Eq. (j3.39p . we obtain 

A., - = -^2Gp{n, - n,) - 2^2eGFq ■ B J ^{h - h) • (3.52) 

The first term on the right side comes just from the background density of matter, and 
the second term is the magnetic field dependent correction. This quantity has been calcu- 
lated for various combinations of temperature and chemical potential of the background 
electrons jSTj IHSl IM| • 

If the denominator of the right side of Eq. ()3.51|) becomes zero for some value of 
Ay^ — Ay^, the value of tan 2^ will become infinite. This is the resonant level crossing 
condition. This was first discussed in the context of neutrino oscillation in a matter 
background by Mikheev and Smirnov 36j, where a particular value of density would 
ensure resonance. Presence of a magnetic field will modify this resonant density, as seen 
from Eq. ()3.52|) . The modification will be direction dependent because of the factor q ■ 3. 
Some early authors [^31111 contemplated that, for large the magnetic term might even 
drive the resonance. However, later it was shown [HZ] that the magnetic correction would 
always be smaller than the other term. So, if one considers values of B which are so large 
that the last term in Eq. ()3.52j) is larger than the first term on the right hand side, it 
means that one must take higher order corrections in B into account. 

It should be noted that the type of corrections to the dispersion relation discussed 
in Sec. 13.41 appear from chiral neutrinos. Thus, they produce chirality-preserving mod- 
ifications to neutrino oscillations. In addition, if the neutrino has a magnetic moment, 
there will be chirality-fiipping modifications as well. Many of these modifications were 
analyzed in the context of the solar neutrino problem, and we do not discuss them here.^ 
As pointed out in chapter ^ they are not important for solar neutrinos, although may 
be important in other stellar objects like the neutron star where the magnetic fields are 
much larger. 

recent paper on chirality-fiipping oscillations is |68| . where one can obtain references to earlier 
literature. Some early references are also found in ^ and j69| . 
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An interesting possible consequence of neutrino oscillations have been discussed (TO] 
in the context of high velocities of neutron stars. In chapter |21 the puzzle about the 
high velocities of the pulsars was addressed in terms of the asymmetric cross-section of 
the inverse beta decay process which compels asymmetric emission of neutrinos from 
the neutron stars. In the present situation the asymmetric emission of the neutrinos is 
generated by resonant conversion of the neutrinos inside the neutron star's body. In a 
material background containing electrons but not any other charged leptons, the cross- 
section of i^e's is greater than that of any other flavor of neutrino. If z/g's can oscillate 
resonantly to any other flavor, they can escape more easily from a star. In a proto-neutron 
star, the resonant density at an angle 9 with the magnetic field occurs at a distance 
i?o + 5 cos 9 from the center, where 5 is a function of the magnetic field and specifies 
the deformation from a spherical surface. So this distance is direction dependent, as 
we have seen from the previous discussion on the resonant level crossing condition in a 
magnetized medium. Therefore the escape of neutrinos is also direction dependent, and 
the momentum carried away by them is not isotropic. The star would get a kick in the 
direction opposite to the net momentum of escaped neutrinos. This was suggested by 
Kusenko and Segre ^Uj , who estimated that the momentum imbalance is proportional to 
5 and can have a magnitude of around 1% for reasonable values of B. Later authors [7T] 
criticized their analysis and argued that the effect was overestimated by them, because 
the kick momentum vanishes in the lowest order in 5. A recent and detailed study [72] 
indicates that these criticisms may not be well-placed, and the kick momentum might 
indeed be proportional to the surface deformation parameter 5. 

3.6 Conclusion 

The main points discussed in this chapter are related to the propagator of a charged 
particle and the neutrino self-energy in a magnetic field and in a magnetized medium. The 
Schwinger propagator acts like the central point related to which all the discussion on this 
chapter follows. This chapter also acts as a precursor for the rest of the material in this 
thesis because what follows heavily relies on the expression of the Schwinger propagator 
in a magnetized medium. 

While discussing about the Schwinger propagator due attention has been paid to the 
phase factor accompanying it. It has been pointed out that the origin and the purpose of 
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the phase factor is related to the gauge degree of freedom of the background gauge field. 
Subsequently the methods of statistical field theory has been utilized to write down the 
expression of the Schwinger propagator in a thermal medium. 

As a magnetic field interacts with the intermediate virtual charged particles, so the 
neutrino self-energy also gets modified from its vacuum value in presence of a magnetic 
field. In the section on neutrino self-energy the expression of the self-energy is derived 
from a purely general idea based on a form-factor analysis. From the discussion on the 
neutrino self-energy it is evident that the general form-factor analysis in a magnetized 
medium yields considerable information. The power of such an analysis is felt when it 
predicts the dependence of the self-energy on the angle between neutrino propagation 
direction and the external magnetic field in a magnetized medium without the help of 
any actual loop calculation. 

Section 13.51 discusses about neutrino oscillations in a magnetized medium. The im- 
portant question about the effect of an external magnetic field on the resonant conversion 
of neutrinos in a medium has been addressed, which is a very important topic in astro- 
physics. In this section another attempt has been made to explain the puzzle about the 
high velocities of the proto-neutron stars by the mechanism of resonant conversion of 
neutrinos from one flavour to the other. In literature this type of mechanism by which 
the proto-neutron stars get a kick in one specific direction due to asymmetric emission of 
neutrinos is often called the 'neutrino rocket mechanism'. 
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CHAPTER 4 



Electromagnetic interactions in a 
magnetized medium 



4.1 Introduction 

This chapter focuses on actual calculations done using the Schwinger propagator. The 
following section will focus on the issue of the effective electromagnetic vertex of the 
neutrinos in presence of a magnetized medium. But before going into the techniques 
of calculation a brief introduction to the topic of electromagnetic vertex of neutrinos is 
presented here. 

Previously a lot of work has been done concerning the neutrino-neutrino-photon vertex 
in the presence of a background magnetic field [SI] . The vertex arises from the diagrams 
of Fig. 11.11 which contain internal M^-lines. In addition, there is a diagram mediated 
by the Z-boson, as shown in Fig. 14.11 For phenomenological purposes, we require the 
electromagnetic vertex of neutrinos only in the leading order in Fermi constant. It should 
be realized that in this order, the diagram of Fig. II. lb does not contribute at all, since it 
has two W^-propagators. The remaining diagrams, shown in Fig. II. lb and Fig. 14.11 can 
both be represented in the form shown in Fig. 14.21 where an effective 4-fermi vertex has 
been used. The effective 4-fermi interaction can be written as, 

^eff = -y2GF['0(^)7^^^H] [V'(£)7a(^v + ^A75)V'w] , (4.1) 

where L = |(1— 75) is the left handed projection operator. If the neutrino and the charged 
lepton belong to different generations of fermions, this effective Lagrangian contains only 
the neutral current interactions, and in that case gv and qa are equivalent to cy and —ca 
defined in Eq. ()3.44|) . On the other hand, if both u and i belong to the same generation, 
we should add the charged current contribution as well, and use 

gy = Cv + l, gA = -{CA + 1) . (4.2) 
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Figure 4.1: Z-photon mixing diagram contributing to the neutrino electromagnetic vertex. 




Figure 4.2: The neutrino electromagnetic vertex in the leading order in the Fermi constant. 



Many processes involving neutrinos and photons have been calculated using the 4-fermi 
Lagrangian of Eq. ()4.H1 . The calculations simplify in this limit for various reasons. First, 
we do not have to use the momentum dependence of the gauge boson propagators. Second, 
since two charged lepton lines form a loop in Fig. 14.21 the phase factors of the propagators, 
as discussed in appendix IB. 11 cancel each other. 

The background magnetic field, can give rise to many physical processes which are 
forbidden in the vacuum. One such process is the decay of a photon into a neutrino- 
antineutrino pair: 



7 



z/ + z/ . 



(4.3) 



This was calculated using the Schwinger propagator in some very early papers [731 12^1 • 
Assuming two generations of fermions, the decay rate was found to be 



487r3a; 



(4.4) 
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where e^, and uj are the polarization vector, the momentum 4- vector and the energy of 
the initial photon, and the quantity ^£ was evaluated in various limits by these authors. 
F^i, is the dual of the electromagnetic field strength tensor. For example, if a; <^ m^, they 
found 

f^f^V for eB<^mj, 
^, = ^sin'ex \''y\J (4.5) 

U(f) fore^»mi, 

where 6 is the angle between the photon momentum and the magnetic field. No matter 
which neutrino pair the photon decays to, both charged leptons appear in the decay rate 
because of the loop in Fig. 14.11 In the calculation only the axial couplings of the charged 
leptons contribute to the amplitude. 

A related process is Cherenkov radiation from neutrinos: 

i/^z/ + 7. (4.6) 

Again, this is a process forbidden in the vacuum. But a background magnetic field 
modifies the photon dispersion relation, and so this process becomes feasible. The rate of 
this process has been calculated by many authors [731 EH I2H1 ES] , and all of them do not 
get the same result. According to Ref. the rate for the process is given by 

r = ^{g'v + gl){eB)'nsm' 9F{n' sin' 9/eB) , (4.7) 

where fl is the initial neutrino energy and 6 is the angle between the initial neutrino 
momentum and the background field. For large magnetic fields satisfying the condition 
eB ^ a;^ sin^ 9, the function F is given by 

Fix) = 1 \ \ \ . (4.8) 

^ ^ 2 3 24 60 ^ ' 

The modification of this process in the presence of background matter has also been 
calculated [7U] . 

Another process that has been discussed is the radiative neutrino decay 

+ (4.9) 



Unlike the previous processes, this can occur in the vacuum as well when the neutrinos 
have mass and mixing. However, a background magnetic field adds new contributions to 
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the amplitude, and the rate can be enhanced. Gvozdev, Mikheev and Vasilevskaya [77j 
calculated the rate of this decay in a variety of situations depending on the field strength 
and the energy of the initial neutrino. For a strong magnetic field {B ^ Bg), they found 
the decay rate of an ultra-relativistic neutrino of energy Q to be 



where U is the leptonic mixing matrix, 6 is the angle between the magnetic field and the 
neutrino momentum, m is the electronic mass, and 



The curious feature of this result is that this is independent of the initial and the final 
neutrino masses. 

The form for the 4-fermi interaction in Eq. 1)4.11) suggests that the neutrino electro- 
magnetic vertex function can be written as |i78j : 



Here, the term n^^, is exactly the expression for the vacuum polarization of the photon, 
and appears from the vector interaction in the effective Lagrangian. The other term, 
n^^, differs from II^i, in that it contains an axial coupling from the effective Lagrangian. 
This equality is valid even when one has a magnetic field and a material medium as the 
background. Thus, the calculation of the photon self-energy in a background magnetic 
field in matter can give us information about the neutrino electromagnetic vertex in the 
same situation. 

Here we concentrate on calculations done in presence of a medium with a uniform 
background magnetic field. We consider various cases where the charged particles af- 
fected by the magnetic fields remain as virtual particles in the Feynman diagrams of the 
respective processes. The techniques of the calculation involve the use of the propagators 
introduced in chapter El This chapter discusses the electromagnetic vertex function of 
neutrinos and consequently the structure of the two second rank tensors U^u and H^^ 
which are necessary to calculate it, as is evident from Eq. ()4.12|) . H^^, has been calculated 
in a magnetic field and in a magnetized medium by various authors |^ 1211 EOl- In the 
following section only the relevant portions of the results are presented briefly. The main 




(4.10) 




(4.11) 
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discussion is around the calculation of 11^^ in presence of a magnetized medium. For the 
sake of comparison the results of the calculations has been given for various backgrounds 
as vacuum, medium, magnetic field and obviously a magnetized medium. Some comments 
on neutrino-photon scattering in a background magnetic field, which is highly suppressed 
in vacuum, is provided in subsection 14.2.11 Section 14. HI discusses the concept of the ef- 
fective electric charge of a neutrino from a quantum field theoretical point of view. As 
neutrinos do not couple to photons in the tree level in the standard model so naturally 
if the neutrinos acquire some electric charge it must be an effective one. This effective 
charge of neutrinos in a magnetized medium has been calculated to odd orders in the 
background magnetic field in section 14.31 Section 14.41 describes the methods for handling 
ultraviolet divergences in calculations involving background magnetic fields. This chapter 
concludes with a general discussion on the various ideas exposed in the following sections. 

4.2 Neutrino-photon scattering and the electromagnetic vertex 
of neutrinos 

In this section we will discuss effective neutrino photon interactions in a magnetized 
medium [221 and the constituents of the neutrino photon electromagnetic vertex. 
But before taking up the issue of the uwy vertex some comments on the neutrino-photon 
scattering in a background magnetic field follows. 

4.2.1 Some comments on neutrino-photon scattering 

The cross-section of neutrino-photon scattering is highly suppressed in the standard 
model due to Yang's theorem jTHj, which makes the scattering cross-section vanish to the 
first order of the Fermi coupling Gp- But in presence of a magnetic field, neutrino-photon 
scattering can occur and to orders of Gp the cross-section has been calculated [20] • In the 
following paragraphs of this subsection a brief overview on neutrino photon scattering in 
a magnetic field background is supplied. 

Gell-Mann showed jHJ that the amplitude of the reaction 

7 + z/^7-Fz/ (4.13) 

is exactly zero to order Gp because by Yang's theorem [73122] two photons cannot couple 
to a J = 1 state. In the standard model, therefore, amplitude of the above process appears 
only at the level of 1/M^ and as a result the cross-section is exceedingly small |53] . 
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Figure 4.3: 



(a) (b) 
The 1-loop effective vertex for two neutrinos and three photons. 



But there is no such restriction on the couphng of three photons with neutrinos as, 

7 + z/^7 + 7 + i/. (4.14) 

The cross-section of the above process can be calculated from the effective Lagrangian 
proposed by Dicus and Repko jHl]. The diagrams for the two neutrino three photon 
interaction are shown in Fig. 14.31 where Fig. 14.3b shows the contribution from the W 
exchange diagram and Fig. 14.3b shows the contribution from Z exchange. Denoting the 
photon field tensor as F and the neutrino fields by '^{v)^ and integrating out the particles 
in the loop the effective Lagrangian comes out as 

^c. = ^^^^ [5(iV,.F-)(F,,F^0 - 14iV,.r^F.,F-] , (4.15) 
where cy was defined in Eq. fj3.44|) . m is the electronic mass, and 

A^M^ = 9^(V'(,.)7i^^^H) - <9,.(V'h7/.^V'm) ■ (4.16) 

For energies much smaller than the electron mass, this can be used as an effective La- 
grangian to calculate various processes involving photons and neutrinos in the presence 
of a background magnetic field Fy,y. For this, we simply have to write 

= /M. + i^M-, (4.17) 

where now /^j^ is the dynamical photon field, and look for the terms involving F^j,. For ex- 
ample, Shaisultanov calculated the rate of 77 — z/z7 in a background field. Eq. ()4.15|) 
shows that in the lowest order, the amplitude for involving z^e's would be proportional to 

^fB B 
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where Be is the value of the magnetic field defined in Eq. p.7|) . Since the amplitude 
without any magnetic field jH3] is of order it follows that the background field 

increases the amplitude by a factor of order {M]y/mYB/Be, or the rate by a factor 
{Mw/'>TLeY{B/Be)'^. Later calculations ^85^ have extended these results by including other 
processes obtained by crossing, like uV — > 77 and z/7 — > z/7. To obtain higher B terms 
in these cross-sections, one needs the effective Lagrangian containing higher order terms 
in the electromagnetic field strength. Such an effective Lagrangian has been derived by 
Gies and Shaisultanov [8Bj . 

Alternatively, the amplitudes can be calculated using the Schwinger propagator for 
charged leptons. Such calculations for 77 —>■ vV were done some time ago [HZIEni- One of 
the important features of this calculation is that in the 4-fermi limit, the diagram contains 
three electron propagators. In such situations, the phase factor '^{x,x') appearing in the 
Schwinger propagator of Eq. (j3.2p cannot be disregarded^. In the calculation, only the 
linear term in B was retained in the amplitude so that the results are valid only for small 
magnetic fields. However, since no effective Lagrangian was used, the results are valid 
even when the energies of the neutrinos and/or the photons are comparable to, or greater 
than, the electron mass. Later authors [SH] reported some mistakes in this calculation 
and corrected them. 

The previous discussions about neutrino-photon scattering relied on effective La- 
grangians and the use of Schwinger propagator to calculate a 3-point function, but aside 
these scattering phenomena there are a lot of phenomena, discussed in the introduction of 
this chapter, which involves the vv^ vertex in a background magnetic field. To calculate 
the rates of these processes one has to calculate the effective electromagnetic vertex of 
neutrinos in a background magnetic field and in a magnetized medium. This section is 
primarily concerned about the electromagnetic vertex of the neutrinos in a magnetized 
medium. 

For simplicity we consider the background temperature and neutrino momenta to 
be small compared to the masses of the W and Z bosons. We can, therefore, neglect the 
momentum dependence in the W and Z propagators, which is justified if we are performing 
a calculation to the leading order in the Fermi constant, Gp- In this limit the 4-fermi 
interaction is given by Eq. ()4.1|) written in the beginning of this chapter. Moreover if we 

detailed discussion on the phase factor of the Schwinger propagator is presented in subsection l3.2.2l 
of chapter O and appendix IB. II 
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restrict the temperature of the system such that muons and taons are not produced in the 
medium then in the leptonic part of the Lagrangian in Eq. ()4.1|) , '?/'(^) stands for electrons 
alone. For electron neutrinos, 

1 
2 



9y = 7: + 2 sin^ ^^w, 



1 

9a = 



2 



For muon and tau neutrinos, 



9y = + 2sin2 6'w, 
1 

With the interaction Lagrangian as given in Eq. 1)4.11) the uw-f vertex, as shown in Fig. 14.21 
can be written in terms of two tensors. The vector- vector amplitude Iln,y{k) and the 
axialvector- vector amplitude 11^^ (/c). This fact is symbolically written down in Eq. ()4.12|) . 
The vector- vector amplitude term Il^^{k) turns out to be: 

^^,.(A;) = (-ze)^(-l) / ^ Tr b,zS{ph.zS{p')] . (4.19) 

Here and henceforth p' = p + k. From the above equation it can be seen that n^j,(fc) is 
the photon vacuum polarization tensor. The axialvector- vector amplitude 11^^^ (A;) comes 
out as: 

K^ik) = {-^en-l) b,l^^S{p)^^tS{p')] . (4.20) 

Both tensors are obtained by calculating the Feynman diagram given in Fig. 14.21 The 
following subsections discusses about the second rank tensors Ilfj,u{k) and n^y(/c). 

4.2.2 The photon vacuum polarization tensor Ii^y{k) in different backgrounds 

The calculation of the photon vacuum polarization tensor in a magnetized medium 
has been done by previous authors [23 IHO] • In this subsection only a brief summary of 
the previously known results is reproduced. The form of Tl^y{k) to one-loop is given in 
Eq. ()4.19|) and it is also pointed out that 11^^ (/c) is a necessary quantity for understanding 
the neutrino electromagnetic vertex. 

The form of 11^;^ (fc) is restricted by two symmetries. Owing to gauge invariance, 'n.^y{k) 
satisfies the conditions 

Fn^,(A;) = 0, rn^,(A;) = 0. (4.21) 
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In addition, Bose symmetry implies 



(4.22) 



These conditions restrict the form of the possible tensor structure of the vacuum polariza- 
tion tensor. The general form of Ilfj,u{k) in various backgrounds is dictated by Eq. ()4.21|) 
and Eq. ()4.22|) as discussed below. 

In vacuum 

In the vacuum, the tensor 11^,^ (fc) depends only on the momentum vector k. Thus, the 
most general form for n^j,(A;) is given as. 



This form satisfies both the conditions given in Eq. ()4.2H) and Eq. ()4.22|) . 
In a medium 

In a medium the expression of n^,y(fc) differs from that in the vacuum. Although 
n^^(fc) still has to satisfy Eq. KTTl and Eq. (jO^ . the form given in Eq. does 
not follow. This is because Il^^{k) can now depend, apart from the momentum vector 
k^, on various vectors or tensors which characterize the background medium. Even for a 
homogeneous and isotropic medium, there is an extra vector in the form of the velocity of 
its center of mass, m^. In general in place of k"^ the Lorentz scalars will now depend upon 



The above expressions of u and K reduce to the energy and momentum of the photon in 
the rest frame of the thermal medium where is given as in Eq. ()3.18|) . 

In presence of a medium without any magnetic field the most general expression of 
the polarization tensor can be written in terms of form factors along with the new tensors 
constructed out of and the ones we already had in absence of a medium as [5H] . 



(4.23) 



K = - fc2 . 



uj = k ■ u , 



(4.24) 
(4.25) 



u^u{K, u) = Ut t^. + L^, + Up 



(4.26) 
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p + k = p' 




Figure 4.4: One-loop diagram for vacuum polarization. 



Here 

T^,„ = g^,u - , (4.27) 
L,. = ^ , (4.28) 

P,.u = ^s^^^^k'^u", (4.29) 

with 

u^ = g^pu'', (4.30) 

and 

k k 

9i.p = 9,.p - • (4.31) 

In Eq. ()4.2(ij) all form-factors are functions of u, K and other parameters as chemical po- 
tential and inverse temperature (3. In Eq. fl4.2(j|l . Hp is nonzero only if the background 
medium or the interaction Lagrangian of the photon break both P and CP symme- 
tries 



In presence of a background magnetic field 

The vacuum polarization of the photon in a magnetic field has also been calculated [221 
up to one loop but to all orders in the external field and the expression of Il^y{k) in our 
convention, as discussed in section 11.31 is given by 

^Ak) = [{9,uk' - k,k.)No - {glkl - fcjl A:|l)iV|| - {g^^kl + kik^)N,] . (4.32) 

Here A'^o, Nti and A^_l are the form-factors, which are functions of B, kfi, and k'j_. 
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In presence of a magnetized medium 

The discussion on Il^i,{k) in presence of a magnetized medium is presented with some 
more details than the previous ones. This is because the form of Il^^{k) is important for 
the expression of the effective charge of the neutrinos in a magnetized medium. It also 
sets the stage for the calculation of n^^(fc) in a magnetized medium. 

At the 1-Ioop level, the vacuum polarization tensor arises from the diagram in Fig. 14.41 
To evaluate this diagram, one needs to use the electron propagator within a thermal 
medium in the presence of a background electromagnetic field. Rather than working with 
the complicated expression for a general background field, we will specialize to the case of 
a uniform magnetic field. Once this is assumed, the field can be taken in the 2;-direction 
without any further loss of generality. 

The amplitude of the 1-loop diagram of Fig. 14. 41 is given by the expression in Eq. (j4.19|) . 
Only in this case the propagators appearing in Eq. (|4.19|) like S{p) and S{p') must be 
replaced by the appropriate ones suitable in a magnetized medium, Sb{p) and Sb{p') as 
introduced in chapter El Eq. ()3.22|) . As the loop is a 2-point function the phase factor 
drops out of the calculation, as is discussed in appendix IB. 11 From Eq. ()3.22p . we see 
that there are two terms in the propagators Sb{p) in a magnetized medium, the vacuum 
part and the other part which involves background matter distribution. In the trace 
appearing in the right hand side of Eq. ()4.19p there appears a product of two propagators 
Sb{p) and Sb{p') with other quantities in between. This product will give rise to four 
quantities of which one term will contain two factors of S^s from the two propagators. 
It has no importance as far as the background effects are concerned. There will be one 
term containing two factors of S^s, which will contribute only to the absorptive part of 
the vacuum polarization. There will be two terms in the trace each containing a factor 
of and a factor of S^. These terms contribute to the present calculation, as they are 
important for calculating the effective charge of the neutrinos in a magnetized medium. 
The relevant part of Ilf^,y{k) is therefore made up of the two terms inside the trace each 
containing one vacuum part and one thermal part of the total electron propagator. With 
this introduction the expression of n^i,(fc) is written as 




(4.33) 
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Using the form of the propagators from Eq. ()3.3p and Eq. ()3.23|) and manipulating under 
the integral signs the 1-1-component^ of the polarization tensor to odd powers in the 
external magnetic field can be written as [HUj . 



. (4.34) 



t: / ^ n , ^ tan(ei3s) tan(ei3s') , ,, 

tan eSs + p'^n tan eSs' - / , ' ' {p + p] 

tan(eo(s + s )) 

The above form of the vacuum polarization tensor is valid in the rest frame of the back- 
ground medium. ?7_(po) contains the information about the distribution of electrons and 
positrons in the background and is given by 

V-{Po) = Vf{po) - Vf{-Po) , (4.35) 

where ?7f(po) is specified in Eq. fj3.17p . In the above equation the argument of the rip 
function is not p ■ u but only the energy component po? because in the rest frame of the 
medium p ■ u = po due to the form of u in the rest frame, given in Eq. ()3.18p . From the 
form of n^j,(fc) as given above it is seen that it is antisymmetric in its tensor indices. 

4.2.3 The axialvector- vector amplitude 11^^^ (A;) in different backgrounds 

The form of the axialvector- vector amplitude 11^^ (fc) to one-loop is supplied in 
Eq. ()4.20j) . In this regard it must be noted that the importance of this second rank 
tensor arises from the fact that it is an important constituent of the neutrino electromag- 
netic vertex as given in Eq. ()4.12j) . The neutrino electromagnetic vertex to one-loop is 
depicted in Fig. 14.21 While discussing about Il^^^{k) it should be remembered that for the 
electromagnetic vertex i.e., the vertex in the Feynman diagram where the external photon 
couples to the fermion loop in Fig. 14. 2t we have the gauge invariance condition, 

nl{k) = 0. (4.36) 

In the following paragraphs the form of 11^^, (A;) in various backgrounds is discussed. 

^The l-l-component of the propagator is discussed in appendix lO where quantum statistical field 
theory is discussed. 
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In vacuum 

n|^^(fc) violates parity due to the presence of 75 in its expression as given in Eq. (j4.2U|) . 
In accordance with this fact the available vectors and tensors at hand are the following, 

k^, g^y and e^^yXu- (4.37) 

The two point axialvector- vector amplitude 11^^ (fc) can be expanded in a basis constructed 
out of the above tensors. Given the axial-tensor like parity structure of 11^^ (fc) the relevant 
tensors at hand are e^^x^k^k'^ and e^u\„g^'^ both of which are zero. So in vacuum I\^^^{k) 
vanishes. 

In a medium 

On the other hand, in a medium in absence of any magnetic field, we have an additional 
vector u^, i.e the velocity of the centre of mass of the medium. Therefore the axialvector- 
vector amplitude can be expanded in terms of form-factors along with the new tensors 
constructed out of and the ones we already had in absence of a medium. A second rank 
tensor constructed out of them can be e^^apu'^k^ , [211 E] which would verify the current 
conservation condition for the two point function. 

In presence of a background magnetic field 

In this case neither C or P but CP is a symmetry, first we look at the CP transforma- 
tion properties of the axialvector- vector amplitude^. As the — Ay mixing Lagrangian 
is of the form, 

^z-A oc Z^'{k)Yily{k)A''{-k) , (4.38) 

we can specify the CP transformation property of n^^,(/c) knowing the CP transformation 
properties of the fields and A^. Under a CP transformation the various components 
of the tensor transforms as. 



CP : 


-'■■'^00 


-^-^00 ' 




CP : 


: n^. - 




(4.39) 


CP : 


'-hi 


^ '-hi ■ 





^This is a preliminary attempt to get some idea about the tensorial basis of n|^^(fc). Strictly speaking 
this attempt to use CP transformation properties of the axialvector- vector amplitude must be incomplete 
because we know that the CP symmetry is also broken in nature. 
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At first, out of all the tensors available in the present situation, a set of mutually or- 
thogonal 4-vectors are constructed with specific CP transformation properties. In the 
next step these different 4-vectors are combined in pairs to produce a second rank tensor 
which has similar CP transformation properties to that of n|^^(A;). As the 4-vectors whose 
combinations give the second rank tensors were mutually orthogonal to start with, the 
second rank tensors are also mutually orthogonal. The basic task in this section is to list 
all the second rank tensors which have similar CP transformation properties to that of 

Before starting the classification of the various second rank tensors at hand which have 
similar CP transformation properties as that of 11^^ (/c), some preliminary remarks about 
the form-factor analysis is in order. In the beginning of this section we have seen that 
Ii^^^{k) is zero in the vacuum. As a consequence it can be inferred that Ii^^^{k) calculated 
in the presence of a background uniform magnetic field vanishes as the magnetic field goes 
to zero. In a uniform background magnetic field, the vectors and tensors at hand are 

F^,, F^., fcf, k>i. (4.40) 

Ffj^y is the field strength tensor for the uniform background magnetic field and F'^^ is 
its dual given as F^^^ = ^^e^''^'^ Fp^. For a magnetic field directed along the 2;-axis the 
expressions of F12, -F21, F03, F30 are given in Eq. in chapter IHl All other components 

of F^j/ zero. The sign of the 4-dimensional totally antisymmetric tensor is specified in 
Eq. ()3.32j) . Now due to Lorentz invariance the form-factors multiplying the basis tensors 
can be functions of /c^, k\ and even functions of the background magnetic field as 
F^^F'^'^, Fp^yF'^'^FpfjF'^^ and other higher powers. Terms containing an odd number of 
FpyS such as, 

Fx^F'^'F^^ (4.41) 

will yield zero because only F12 and F21 are non zero. The first thing to notice from 
the above discussion is that the form-factors themselves are CP even and so the CP 
transformation property of 11^^ (A;) is solely dependent on the basis tensors of 11^^^ (A;). 
Secondly, from Eq. ()3.3p of chapter El it is evident that in the B ^ limit the propagator 
in a magnetic field transforms into the vacuum propagator and it is shown ^55j that 
the Schwinger propagator can be expanded in a series of powers of the external field 
magnitude i.e.B. From this observation it is assumed that calculations based on the 
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Schwinger propagator are perturbative in the the external field i.e. can be expanded in 
a series of B. In the present circumstance the form-factors can also be expanded in a 
power series of the the even powers of B where the first term of the series is B^ raised to 
the power zero i.e. a term independent of B. So if the external magnetic field vanishes 
the form-factors themselves in general do not vanish. As a result the second rank tensors 
which multiply the form-factors must tend to zero as the external field is put off because 
Il^^^{k) itself vanishes when the external field is absent. 

Now all the possible 4-vectors that can be made from the quantities enlisted in 
Eq. ()4.4()|1 are given. They are, 



b1 = (Fk) 



The expressions (Fk)^, and (Fk)^ stands for 

m, = F,.k^, 

(Fk), = F^X . ^^-^"^ 

The set of 4-vectors 6f s, i = 1, 2, 3, 4 are mutually orthogonal to each other and can serve 
as the basis vectors to build up the tensor basis of li^^^ik). 

Next the CP transformation properties of these vectors are summarized. 

CP : b1 ^b1, . . 

CP : 6i ^ 6i. ^ ^ ^ 

The other three vectors have similar transformation properties as 

CP ■ h Z h ^^-^^^ 

From Eq. (|4.39|) we can see that a suitable tensor basis can be built up from vectors 6f 
where i = 2,3, 4. The CP transformation of the axialvector-vector amplitude compels us 
to disregard bi basis vector. 

Now we can list all the possible second rank tensors, having similar CP transformation 
properties to that of 11^^ (A;), made by combining in pairs the mutually orthogonal 4-vectors 
62, &3, &4- There are nine of them. For later usage we denote them as follows: 

B't=b'ibl={FkYkl, B^''=b^b^^ = {Fkyk['^ , B^''=bi^bl={FkYk'i, (4.46) 

B!;''=b^b^^ = {Fkyk'i , B^''=b^bl=ki^^kl , Bl^'' =b^^bi=kie^ . 
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This basis gives nine second rank mutually orthogonal tensors. The orthogonal condition 
for the second rank tensors is given by, 

fir5iM- = 0. (4.47) 

In the above equation i,j = 1, .., 9 and i ^ j. Any second rank tensor containing higher 
field dependence can be represented by suitable combinations of these tensors. 

Out of these nine basis tensors, Bl^'^ , B^'^ , Bg'^ and Bg'^ do not vanish in the B — > 
limit and are therefore redundant. The remaining five second rank tensors qualify 
successfully as probable building blocks of the axialvector-vector amplitude in the light 
of CP transformation properties. But till now the gauge invariance condition as given in 
Eq. fj4.36|) has not been used. Regarding gauge invariance the first thing to note is that 
Bi" can also be written as, 

= {FkYiFky 

= -B\gfkl - A;f A;p , (4.48) 

and as a result 

A;^5f = 5f A;^ = . (4.49) 

From the above equation it is seen that Bi'^ satisfies the gauge invariance condition for the 
photon vacuum polarization tensor. From Eq. ()4.48|) and the expression of the vacuum 
polarization given in Eq. ()4.32j) it can be seen that, except a numerical constant which 
can be absorbed in N\\, B^^ actually occurs in the expression of the polarization tensor. 
From these observations it can be inferred that B^^ is not a suitable tensor to be in 
the set of the probable second rank tensors which can act as the basis of axialvector- 
vector amplitude. The basis tensors of the axialvector-vector amplitude can be made by 
suitable combinations of the other four second rank tensors listed above, satisfying the 
CP transformation properties as given in Eq. ()4.39|) and the gauge invariance condition. 
They are, 

Pf- = {Fkyk'^^ , (4.50) 

and 

p^"" = [FkYki + {Fkye^ + -A[{FkYk\ - [Fkyk^] . (4.51) 
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P2'^ can be simply written as: 

= {FkYkl + {Fkyk^l + klF^"" , (4.52) 

utilizing the identity, 

F^"^ = ^ UFkfk^^ - {Fkykff] . (4.53) 
The result as given in the papers by Hari Dass and Raffelt ^ verifies this choice 



=3 



Kik) = f\\Pi,^ + C±P2,.\ , (4.54) 

where Cy and C± are the form-factors. 

In presence of a magnetized medium 

In presence of a magnetized medium the situation complicates. In this analysis we are 
not going into an in depth study of the tensorial basis as is done in the case where there 
was no medium. No new techniques are used in this discussion, it is only a continuation 
of the previous analysis in a new situation. 

To start with we again emphasize on the CP transformation properties of the 
axialvector- vector amplitude 11^^ (A;). Unlike the vacuum case now the theory may not 
be CP invariant. If the background contains more number of electrons than positrons 
and/or if the electrons momenta are not random but along some specific direction then 
CP is violated by the background. For simplicity here we will discuss only those cases 
where the background does not break CP. 

Now the form-factors can be functions of odd powers of the magnetic field as well, as 
now new scalars as: 

iFk)u = {Fk)X , u 
{Fk)u = {Fk)^u^ , ^ ' 

are also available. These scalars change sign under CP transformation. Some of the 
form-factors may contain odd powers of the scalars listed above multiplied by equal odd 
powers of chemical potential of the background charged fermions. These form-factors will 



^However the metric used by the authors in references mentioned is different from that of ours, and in 
their calculation the /i vertex is the vector type vertex. Due to this change the sign of fc^F'^" appearing 
in P!^'^ is different from the one they obtained. The metric used in this thesis is specified in section [1.31 
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then not change sign under a CP transformation. So in a magnetized medium there can 
be two kinds of form-factors, one which does not change sign and the other which does 
change sign under a CP transformation. This is unhke the previous case where there 
was no medium. So in this case the basis tensors can have different CP transformation 
properties as the form-factors which muhiply them can also have different transformation 
properties. The vector bi in Eq. (j4.42j) had a CP transformation property which was 
different from the other three vectors &2 3 45 a result it was dropped from the set of 
vectors which could combine in pairs to give a set of second rank tensors having similar 
CP transformation property as that of Il^^^{k). In the present circumstances bi is not 
excluded as in the vacuum because the CP transformation property of the basis tensors 
have changed. 

In presence of a medium we can have two different sets of mutually orthogonal vectors. 
The first set is as supplied in Eq. (j4.42|) . The other set of orthogonal vectors useful in a 
medium, are 



In this list we have omitted two vectors. One is {FuY and the other one is m^^. The reason 
is that ultimately we are interested in the rest frame of the medium. In the medium rest 
frame there is no electric field. Also -u^ = in this frame. 

This above set of vectors as given in Eq. ()4.56|) has similar CP transformation prop- 
erties to those of 62 and 63. But the vectors given in Eq. ()4.42j) are not orthogonal to the 
vectors specified in Eq. fl4.56p . So although the number of vectors at hand has increased, 
four from Eq. ()4.42j) and two from Eq. ()4.56|) . the number of mutually orthogonal vectors 
have not increased. So this six vectors as such cannot serve as the basis vectors to build 
up the second rank tensor basis of Tl^^^,{k). Only a linear combination of the two sets 
of vectors can make a orthogonal vector basis. Now we list the set of orthogonal basis 
vectors which can be made from the two set of vectors, they are: 



In a magnetized medium we have these four basis vectors which serves as the building 
blocks of the axialvector-vector amplitude. The basis tensors in this case will be the 



b'l" = {FuY , 



(4.56) 



bl^" = [fuY + {FkY , 
bT = kl , 

b'i'' = kif + u':. 



(4.57) 
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direct product of these basis vectors. There will be sixteen of them but all of them will 
not be useful. The electromagnetic current conservation condition will reduce the number 
of admissible basis tensors. 



4.2.4 One loop calculation of the axialvector- vector amplitude 11^^, (A;) in a 
magnetized medium 

Since we investigate the case with a uniform background magnetic field, without any 
loss of generality it can be taken to be in the z-direction. The relevant Feynman diagram 
of the neutrino electromagnetic vertex to one loop appears in Fig. 14.21 Following that 
diagram the axialvector- vector amplitude 11^^ (A;) is given as in Eq. (|4.2LH) . For calculating 
n^^(A;) in presence of a magnetized medium the propagators appearing in Eq. (j4.2(jp like 
S{p) and S{p') must be replaced by the appropriate ones suitable, Sb{p) and Sb{p') as 
introduced in chapter El Eq. ()3.22|) . 

Similar to the case of n^j,(fc) here also the phase factor does not appear because the 
quantity in question is a one loop 2-point function. From Eq. ()3.22|) . we see that there 
are two terms in the propagator Sb{p) in a magnetized medium, the vacuum part and the 
other part which involves background matter distribution. In the trace appearing in the 
right hand side of Eq. ()4.20|) there appears a product of two propagators Sb{p) and Sb{p') 
with other quantities in between. This product will give rise to four quantities of which 
one term will contain two factors of SbS from the two propagators. It has no importance 
as far as the background effects are concerned. There will be one term containing two 
factors of SbS, which will contribute only to the absorptive part of the axialvector-vector 
amplitude. There will be two terms in the trace each containing a factor of Sb and a 
factor of Sb- These terms contribute to the present calculation, as they are important for 
calculating the effective charge of the neutrinos in a magnetized medium. 

With this introduction the relevant expression of n^j^(/c) can be written as, 

^K^ik) = / ^ Tr [l,l,SUph.^Slip') + l,l,^Sliph.Slip')] . (4.58) 

Using the form of the fermion propagators in presence of a magnetic field and a magnetized 
medium, given by Eq. ()3.3|) and Eq. ()3.23|) in sections 13.21 and 13.31 of chapter El we get 
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+ Tr [7m75G(-p', s')-f^G{-p, s)] Vri-p ■ u)] 
= -e'l -^.j^J^^"^'''^ fds'e''^''^'^R,.{p,p',s,s'), 



(4.59) 



where R^iy(p,p', s, s') contains the traces. 

to even and odd orders in magnetic field 

We calculate R^j,(p, p', s, s') to even and odd orders in the external magnetic field 
and call them and Rj,°j. The reason for doing this is that the two contributions 
have different properties as far as their dependence on medium is concerned, and the 
contributions are 



Rg = 4i77_(p • u)e^,^ap (1 + tan(eSs) tan(eBs')) 

+ sec'(eBs') + p^V^" sec'(eBs) 



and 



R^^j — '^ir]+{p ■ u) m^£^,^i2(tan(eBs) + tan(eBs')) 

+ (fi'M«||P""pI.± + 5'i.a||P""p'^±) sec^(ei3s')} tan(ei3s) 

+ {(i/Ma||P'""Pi^|| -«a||P'"" + yi.a||P'""PM|| ) 

+ (5Ma||P'""Pz.± + 5z.a||P'""P/.±) sec^ {eBs) } tan(eBs') 



(4.60) 



(4.61) 



Here 



ri+{p ■ u) = r]F{p ■ u) + r]F{-p ■ u) , 
rj^{p ■ u) = r]F{p ■ u) - rjF{-p ■ u) , 



(4.62) 
(4.63) 



which contain the information about the distribution functions. Also it should be noted 
that, in our convention 
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(4.64) 



As stated we have split the contributions to Il^^^{k) to odd and even orders in the external 
constant magnetic field. The main reason for doing; so is the fact that, Ill^°\k) and 
]l^^\k), the axialvector- vector amplitude to odd and even powers in eB, have different 
dependence on the background matter. Pieces proportional to even powers in B are 
proportional to i]^{p ■ u), an odd function of the chemical potential. On the other hand 
pieces proportional to odd powers in B depend on ri+{p ■ u), and are even in /i and as a 
result it survives in the limit ^ 0. 

From Eq. (j4.6Up we notice that Il^^^{k) to even orders in the magnetic field satisfies 
the current conservation condition in both the vertices. In Eq. ()4.(iHl we see that all the 
terms in the right hand side are symmetric in the fi and u indices except the first term. 
This term differentiates between the two vertices in this case and as 11^^ (/c) to odd orders 
in magnetic field is gauge invariant in the u vertex we do not get the same condition for 
the axialvector vertex. If in Eq. (|4.61|) we put m = then all the terms in the right will 
be symmetric in both the tensor indices, and as a result current conservation condition 
will hold for both vertices. If the mass of the looping fermion is not zero then from the 
above analysis we can say that only Eq. (j4.36j) will hold. 

If we concentrate on the rest frame of the medium, then p-u = Po- Thus, the distribu- 
tion function does not depend on the spatial components of p. In this circumstance we can 
manipulate under the momentum integral signs, and modify the above equations. The 
techniques used to convert expressions inside the momentum integrals is briefly discussed 
in appendix ^ while other related calculations showing the derivations of the equations 
given below are supplied in appendix |F| Using the relations discussed in the appendix we 
can rewrite, 

R^rj = Airi.ipo) [£^.„„^||P°iip"'ii(l + tan(ei3s) tan(ei3s')) + £;..a||/3±p""p'^^ sec^{eBs') 
+ £^.a±/3,|P°V^" sec^eBs)] , (4.65) 



and 

RS = 4277+ (po 



—e 



sec'^(eBs)tan^(eBs') , , , ^ . , ^ 

^ ' ^ ^A;i + (A;-p)||(tan(eSs) +tan(eSs')) 



tan(ei3s) + tan(ei3s' 

+ 2£:^i2a||(p'j,||P°" tan(ei3s) tan(eSs')) + g^j^a^^ky^ (tan(ei3s) - tan(eSs')) 

"iip^_L sec (ei3s) tan(ei3s') . (4.66) 
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The sign '=' means that the expressions on both sides of the above equations, though 
not necessarily equal algebraically, yield the same momentum integral. To be precise if 
there are two functions of 4-momenta p, p' and parameters s, s', as A{p,p' , s, s') and 
B{p,p', s, s') and 

A{p,p',s,s') = B{p,p',s,s') , (4.67) 

then in general, 

A{p, p, s, s') ^ B{p, ps, s') , (4.68) 

but 

J d^pA{p,p', s, s') e*(f'^)+*(p''«') = J d^pB{p,p', s, s') e*(P'^)+*(p''^') . (4.69) 

As the momentum integrals appearing in the right hand side of Eq. (|4.59|) are of the 
same form as those in Eq. ()4.69|) we can safely use the right hand sides of Eq. ()4.65j) and 
Eq. (|4.66p as representatives of R^fj and R^°j . 

Before going into the next section we comment on the nature of the integral appearing 
in Eq. (j4.59|) . The first point to make is that from the form of R|fj in Eq. (|4.65|) we note 
the axialvector-vector amplitude in a magnetized medium to even orders in the magnetic 
field is antisymmetric, as it was in a medium without any magnetic field. Contrary to 
this R|j°j does not have any well defined symmetry property. 

Secondly, as the integrals are not done explicitly something must be said about possible 
ultraviolet divergence which may creep up in evaluating them. It is a complicated topic 
but fortunately the integrals we are dealing with do not have any divergences as we 
are working in presence of a medium which offers a natural ultraviolet cut-off as the 
temperature of the system. The penultimate section in this chapter discusses the points 
about ultraviolet divergences of the two cases we have considered in this chapter. 

The integral expression for Il^^^{k) in our case does not have any singularities. We can 
now write the full expression of the axialvector-vector amplitude as 

where Rjfj and Rg are given by Eqs. ()4.66p and ()4.65p in the rest frame of the medium. 
A thorough proof of the gauge invariance of n^j,(A;) is given in appendix [Gl 
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4.3 Neutrino effective charge in various backgrounds 

4.3.1 Definition of eff'ective charge 

Intuitively when a neutrino moves inside a thermal medium composed of electrons and 
positrons, they interact with these background particles. The background electrons and 
positrons themselves have interaction with the electromagnetic fields, and this fact gives 
rise to an effective coupling of the neutrinos to the photons. Under these circumstances 
the neutrinos may acquire an "effective electric charge" through which they interact with 
the ambient plasma. The proper definition of this charge follows. 

The off-shell electromagnetic vertex function Fj, is defined in such a way that, for 
on-shell neutrinos, the uwy amplitude is given by: 

M = -iu{q')rM(l)A"ik), (4.71) 

where, k is the photon momentum. Here, u{q) is the the neutrino spinor and A'^ stands 
for the electromagnetic vector potential. In general r,y would depend on k and the char- 
acteristics of the medium. With our effective Lagrangian in Eq. ()4.1|) . the form of T^, is as 
given in Eq. ()4.12|) . The effective charge of the neutrinos is defined in terms of the vertex 
function by the following relation [211 E] : 

Ces = ^ u{q) ro(fco = 0, fc ^ 0) u{q) . (4.72) 

For left handed massless Weyl spinors this definition can be rendered into the form: 

eeff = 7^ Tr [Fo(fco = 0, fc ^ 0) (1 - 7') A • (4.73) 
2go 

Both the vector-vector amplitude and the axialvector-vector amplitude contribute to 
the effective charge of neutrinos. In the following subsections the the effective charge of 
neutrinos in different backgrounds are discussed. 

4.3.2 In vacuum 

In vacuum 11^,^ (A;) vanishes in the limit /cq = 0, — > and 11^^ (A;) is zero. As a result 
there can be no effective charge of neutrinos in vacuum. 

4.3.3 In a medium 

In a medium the axialvector-vector amplitude is of the form e f^^aisu'^ k^ , and so it does 
not contribute for the effective electric charge of the neutrinos since for charge calculation 
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we have to put the index z/ = 0. In the rest frame only the time component of the 4- 
vector u exists, that forces the totally antisymmetric tensor to vanish. In a medium the 
vector-vector amplitude is same as the vacuum polarization as given in Eq. (j4.26|) except 
one point. Eq. ()4.2fi|) contains Up which is an effect of P and CP violation, either by the 
Lagrangian or by the background medium. Phenomenologically it will come in orders of 
Gp and so its inclusion will make the effective charge of neutrinos to be higher in orders 
in the Fermi coupling. So the effect of Hp is neglected in the effective charge calculations 
of neutrinos. The longitudinal projector L^^, is not zero in the limit = 0,k and 
Ul is also not zero in the above mentioned limit [22 El- This fact is responsible for giving 
nonzero contribution to the effective charge of neutrinos in a medium. The effective charge 
of neutrinos in a medium has been calculated previously by many authors [5^ 1^ IHj. 

4.3.4 In presence of a background magnetic field 

From the Eq. (j4.54j) we see that the axialvector-vector amplitude in a background 
magnetic field without any medium does not survive when the momentum of the external 
photon vanishes. Also from Eq. (j4.32|) it is evident that U^^lk) also vanishes when the 
momentum of the external photon vanishes. As a result there cannot be any effective 
electric charge of the neutrinos in a constant background magnetic field. Actually this 
formal statement could have been spoilt by the presence of possible infrared divergence 
in the loop; i.e to say in C\\ and C± p8| and A'o, N\i and A^_l [23]. Since the particle inside 
the loop is massive so there is no scope of having infrared divergences. 

4.3.5 In presence of a magnetized medium 

To odd orders in the external magnetic field the expression of n^j,(A;) is given by 
Eq. ()4.34|) and it shows that the result is dependent on k^. In the limit when the momen- 
tum of the external photon goes to zero 11^;^ (A;) will vanish unless the parametric integrals 
can produce terms which cancels k^. In the present circumstance this cancellation does 
not occur and as a result Ilf^,^{k) vanishes. On the other hand all the components of 
n^j^(A;) does not vanish in the specified momentum limit. The relevant components of 
n^i^n^) which are required for the effective charge calculation, vanish completely when 
the external momentum goes to zero. But all the components of Il^^°\k) does not vanish 
when the external photon momentum goes to zero. As a result this can give rise to an 
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effective neutrino cliarge in a magnetized medium. Tlie calculations showing the validity 
of the above statements appear in appendix |F| 

Now we concentrate on the zero momentum limit of that part of the axial polarization 
tensor which is going to contribute for the neutrino effective charge in the rest frame of a 
magnetized medium. 

Effective charge to odd orders in external field 

Denoting n^^,(fco = 0, fc ^ 0) by 11^^, we obtain 

X r]+{po) [2pl - (A; ■ £^012 (4.74) 

the other terms turns out to be zero in this limit. Although it looks that the term 
[k ■ will vanish in the limit when ko = and k —>■ 0, but it does not. This is 
because the parametric integrals yield terms which contain {k ■ p)\\ in the denominator of 
the momentum integrals and there remains a scope for cancellation. The calculation in 
appendix IF] clarifies the point. The above equation shows that, except the exponential 
functions, the integrand is free of the perpendicular components of momenta. This is a 
peculiarity of this case that the perpendicular excitations of the loop momenta are only 
present in the phase like part of the integrals and in effect decouples from the scene once 
they are integrated out. Its presence is felt only through a linear dependence of the 
external field B when the perpendicular components of k vanish. Upon performing the 
gaussian integration over the perpendicular components and taking the limit ^ 0, we 
obtain. 



lim ^^^F^ / 7^ / dse'''\r ds'e 
^° fco=o,fc^o A-n J (27r)2 i_oo Jo 

X r7+(po) [2p^- (A;-p)||] . (4.75) 

The details of the calculation for the effective charge is given in appendix |Fl For a 
classical gas background the expression of the effective charge, to odd orders in the external 
magnetic field, is given by, 

cos^ cosh(/3/i)-ft'i(m/3) . (4-76) 

Here 9 is the angle between the neutrino three momentum and the background magnetic 
field. The superscript v in e^g denotes that we are calculating the effective charge of 

72 



the neutrino. Ki{mj3) is the modified Bessel function (of the second kind) of order one 
which sharply falls off as we move away from the origin in the positive direction. As 
temperature tends to zero, Ki{m[3) grows as e~™'^, as a result from Eq. (j4.76|) it is seen 
that the effective charge vanishes. In [HSl a more general treatment of the effective charge, 
relinquishing the assumption of a classical gas, has been presented. 

Before ending this section a general discussion can me made on the form-factors of 
n|^j^(A;) in the zero momentum limit. In a background magnetic field the form-factors can 
be functions of scalars containing the magnetic field. The scalars can be of the following 
form: 

k'^F^.F^^k^ and (4.77) 

or 

(FuYiFu)^ and {FuYk^. (4.78) 

The scalars listed above are the most elementary ones, there can be many more possible 
scalars where more number of fc^s, F^'^s and u^s are mutually fully contracted. The 
thing which must be noted is when k tends to zero only terms that can survive in the 
form-factors must be even functions of B. 

Of all possible tensorial structures for the axialvector- vector amplitude in a magnetized 
plasma, there exists one which is independent of the external momentum k given by, 

F^^u^'ul. (4.79) 

It is seen from the above expression, which is odd in the external field, that it survives 
in the zero external momentum limit in the rest frame of the medium. We have earlier 
noted that the form-factors which exist in the rest frame of the medium and in the zero 
momentum limit are even in powers of the external field. This tells us directly that the 
axial polarization tensor must be odd in the external field in the zero external momentum 
limit, a result which is verified by actual calculations. Some similar calculations regarding 
the axialvector- vector coupling is also done by Konar and Das pU] . 

4.4 A short note on divergences in a background magnetic field 

While calculating closed loops in a Feynman diagram in vacuum we encounter ultra- 
violet divergences in many places. For renormalizable quantum field theories there are 
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standard methods of handling this divergences. First the theory is regulated and then 
some of the bare parameters in the theory are rescaled to give a meaningful theory. There 
are various regularization schemes as the cut-off regularization, the Pauli-Villars regular- 
ization and the most prominently used dimensional regularization. Not going into the 
full details of the renormalization programme in vacuum, which has been studied exten- 
sively, here we focus on the methods by which ultraviolet divergences are handled in those 
calculations where we have a background magnetic field. 

The first thing to note here is that magnetic fields do not bring in any new kind 
of divergences, they only modify the properties of the divergences present in the vacuum 
structure of the theory. Vacuum polarization of photon contains ultraviolet divergence and 
consequently the divergence remains when we calculate it in a background magnetic field. 
The expression of Yl^^^^k) calculated in a magnetic field includes the B = vacuum result 
implicitly, which we denote by limg^o n^,y(A;). The divergent part in 11^,^ (fc) calculated in 
a background magnetic field is actually present in lim^^o n^jy(A;). The usual procedure is 
to fix the divergent part in lim^^o ^^iu{k) and subtract it from the the final expression of 
^fMu{k) calculated in presence of a magnetic field, to get rid of the divergence. 

The ultraviolet nature of n|^^(A;) is interesting. In vacuum we have seen from subsection 
14.2.31 that n^j^(A;) is zero and so it cannot have any divergences when we calculate it 
in presence of a background magnetic field. In this connection it can be said that in 
absence of the medium but in presence of the background magnetic field another divergent 
structure could arise, the Adler anomaly [HI], due to the presence of the axial- vector vertex 
in the Feynman diagram of n|^^(A;). A brief discussion of this topic is given in the work 
by Raffelt and loannisian j2H]- 

As stated previously a medium affects the ultraviolet nature of the calculations. The 
finite temperature of the system automatically act as an ultraviolet cut-off of the system 
and so momentum integrals cannot run to infinity to damage the logical consistency of 
the theory. As a result the calculations in a magnetized medium are also expected to be 
free from ultraviolet divergences. 

4.5 Conclusion 

This chapter was about the electromagnetic interactions of neutrinos in a magnetized 
medium. In the introduction various decay processes as the photon splitting into two 
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neutrinos, the Cherenkov process and radiative decays of heavy neutrinos were discussed. 
Section 14.21 deals with neutrino-photon scattering and the electromagnetic vertex of neu- 
trinos. The topic of neutrino-photon scattering, which is highly suppressed in the standard 
model, is briefly discussed in subsection 14.2.11 The interesting thing about this scattering 
process is that in presence of an external magnetic field the neutrino-photon cross-section 
increases which can be an important ingredient for astrophysics. 

In the standard model of particle physics neutrinos do not interact with photons in 
the tree level so the neutrino electromagnetic vertex is an effective one mediated by the 
charged fermions in the loop. Section 14.21 deals with the two important second rank 
tensors, 11^,^ (fc) and n^^^(fc), which are the building blocks of the electromagnetic vertex 
function of the neutrinos. The structure of Il^^{k) in various possible backgrounds has 
been discussed in subsection 14.2.21 Here for the first time methods of statistical field 
theory is used to find out the form of n^j,(A;), to odd orders in the external field, in a 
magnetized medium. Subsection 14.2.31 of the present chapter deals with the axialvector- 
vector amplitude 11^^^ (A;) in a magnetized medium and its possible tensor structures in 
various possible backgrounds. Section 14.31 deals with the issue of the effective charge 
of the neutrino in a magnetized medium. The formula of the effective charge calculated 
shows that it depends upon the angle between the neutrino 3-momentum and the external 
magnetic field directions. Also the effective charge vanishes as temperature tends to zero. 

As no discourse in quantum field theory can be complete without solving the thorny 
issue of divergences popping out from the loops, in section lOl this issue is discussed. It is 
pointed out that the temperature of the medium can act as an ultraviolet regulator for the 
various processes which we have discussed and we can safely assume that the calculations 
are free of ultraviolet divergences. 



75 



CHAPTER 5 



Conclusion 

It is now an established fact that magnetic fields of various magnitudes and various 
shapes pervade our universe. They are present in the intergalactic medium, in the galactic 
region and also in planetary and stellar atmospheres. Order of magnitude estimates of 
these fields have been presented in chapter 

Two different ways have been followed in this thesis to incorporate the magnetic 
field effects. The first way is straight forward and involves the solution of the Dirac 
equation in presence of a classical uniform background magnetic field. In presence of 
a magnetic field the Dirac equation yields an exact solution. The exact wave functions 
in a uniform magnetic field can be used to produce a consistent quantum field theory 
and various scattering cross-sections involving charged fermions can be calculated, as 
elaborately discussed in the chapter |21 Magnetic fields found in astrophysical objects like 
neutron stars, active galactic nuclei vary both in magnitude and direction but the typical 
length scale D in which this variation occurs is much much higher than the Compton 
wavelengths of the elementary particles and as a result this variation can be neglected in 
a quantum field theoretical calculation^". 

Magnetic fields also enter into the elementary particle regime indirectly through the 
effect of virtual charged particles. To evaluate quantities like neutrino self-energy and vac- 
uum polarization of photon intermediate virtual charged particles like electrons, W-bozons 
are required whose propagators are modified in a magnetic field. Although ultimately the 
intermediate particles are integrated out, the magnetic field effects remains frozen in and 
affect the final results. In chapters El and |3] we have encountered some of the calculations 
involving the Schwinger propagator and seen how magnetic fields creep in to modify the 
physical properties of those particles which does not have any interaction with a magnetic 
field normally. 

lOpQj. neutron stars D ^ lOKm. 
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Historically the URCA processes, as specified in section IT^ were the first to be studied 
in an external field. Closely related with the URCA processes is the inverse beta decay 
process. In chapter |21 the inverse beta decay cross-section of arbitrarily polarized neutrons 
in an external magnetic field is calculated. The magnitude of the external magnetic field 
assumed is much less than m^/e and so the magnetic field modification of the proton 
wave function is neglected. Using the solutions of the Dirac equation for electrons in a 
background of a magnetic field a quantum field theoretical formalism is built up to tackle 
various scattering processes consistently. From Eq. (j2.76j) of chapter |21 it is seen that the 
inverse beta decay scattering cross-section in a background magnetic field has a smooth 
i3 ^ limit. In Eq. (j2.7!2j) the cross-section of the inverse beta decay process is given and 
is seen to be dependent on the angle between the incoming neutrino 3-momentum and 
the magnetic field vector. From Fig. 12. H it is apparent that the cross-section of the inverse 
beta decay in an external magnetic field becomes considerably enhanced compared to its 
corresponding value in vacuum when the magnetic field magnitude is greater than Se, 
the critical field defined in Eq. (|1.7|) . Moreover when the initial neutrinos are taken to be 
monochromatic Fig. 12.11 shows that the cross-section of the inverse beta decay process is 
plagued with spikes which goes all the way to infinity. This disease is cured and the spikes 
smear out, as shown in Fig. 12.21 when a flat probability distribution of initial neutrino 
energy is considered. Using the fact that the beta decay cross-section is sensitive to the 
angle between the neutrino 3-momentum and the magnetic field direction attempt has 
been made to explain the high velocities of pulsars, of the order of 450 ± 90 Kms~^ in 
section 12.51 

The charged fermion propagator in a background magnetic field as derived by 
Schwinger is given in chapter 01 subsection 13.2.21 The propagator is defined as an integral 
over a parameter called 'proper time'. Using form-factor analysis the general structure 
of the neutrino self-energy in a magnetized medium has been found out in section 13.41 
From Eq. ()3.39|) it is seen that in a magnetized medium the neutrino self-energy becomes 
sensitive to the angle between the neutrino propagation direction and the magnetic field 
direction. It is also sensitive to the nature of the background i.e., whether the background 
is made up of electrons or muons etc., as is evident from Eq. ()3.42j) . Both of these facts 
affect the resonant level crossing condition in context of neutrino oscillations discussed in 
section 13.51 
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The topic of electromagnetic interaction of neutrinos is discussed in chapter 0] All the 
calculations assumes an effective 4-fermi vertex. There are various processes which are 
either forbidden or highly restrained in vacuum, like a photon decaying into a neutrino- 
antineutrino pair, Cherenkov radiation of neutrinos, neutrino-photon scattering, which 
can become possible in presence of a background magnetic field. These processes, briefly 
discussed in sections 14.11 and 14.21 opens up new channels for energy emission from newly 
born neutron stars and so are astrophysically important. The neutrino electromagnetic 
vertex function can be written in terms of two second rank tensors Il^u{k) and n^^,(fc) 
as mentioned in chapter ^ and explicitly shown in Eq. ()4.12j) . n^j,(fc) is exactly the 
vacuum polarization of the photon and its form in a thermal medium and in presence of a 
magnetized medium is presented in subsection 14.2.21 The general form of the axialvector- 
vector amplitude Il^^^{k), as defined in Eq. (j4.2(J|) . is calculated in a magnetized medium 
in section 14.2.31 A general tensor analysis based on the discrete symmetries shows that 
n^^(A;) must vanish in vacuum. A similar analysis gives the tensor structure of Il^^^{k) in 
a thermal medium. Based upon CP transformation property of n^j^(A;) the tensor basis 
of the axialvector-vector amplitude is analyzed in a background magnetic field. Some 
comments on the tensor basis of 11^^ (fc) in a magnetized medium is supplied at the end 
of subsection 14.2.31 Subsequently Il^^j^{k) has been calculated in a specific limit i.e., when 
ko = and k 0, which is related to the effective charge of the neutrinos. In this limit it 
is observed that 11^^^ (A;) does exist to odd orders in the external magnetic field but to even 
orders it does not. It is seen that when ko = and k ^ 0, n^j,(fc) becomes linear in the 
external field. Using a classical gas approximation the expression of the effective charge to 
odd orders in the external field has been written down. The formula of the effective charge 
calculated shows that it depends upon the angle between the neutrino 3-momentum and 
the external magnetic field directions. Also the effective charge vanishes as temperature 
tends to zero. 

To conclude it must be said that most of the effects discussed in this thesis depends 
upon large magnetic fields, where B > Be- These fields are obtained only at astronomical 
distances, presumably in neutron stars and magnetars. For such distant objects, obser- 
vational data are not clean enough to resolve the effects of the magnetic field. Perhaps 
more accurate measurements of indirect evidences like the velocities of the pulsars or rate 
of cooling of newly born neutron stars can shine some light on the calculations done in 
this thesis in near future. 
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APPENDIX A 



Spin Sum in a uniform background 

magnetic field 



From Eq. ()2.24|1 the spin sum J2sUs{y,n,j\)Us{y*,n,Pi^) can be written as: 

1 



J2Us{y,n,p^)Us{y^,n,j\) 



En + m 



(A.l) 



ij=n~l 

where the Ti j s are 4x4 matrices. The Tj j matrices are obtained by performing the 
matrix multiphcation in the left hand side of the above equation. 

Using the dispersion relation = pi + 171"^ + 2neB, Tn^n can be written as, 



T 



In the 2x2 notation the above matrix can be written as, 
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(A.2) 



T = K 2 
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(A.3) 



where (T3 is the third Pauli matrix. In the 4x4 notation Eq. ()A.3|) can be written as, 

1, 



T 

■J- 7 



n,n - 2 L^(l - ^z) + Enil" + 7^') " Pzil'l" + Y)] 



;[m(l-(Tj+|^|| +|^||75], 



(A.4) 
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(A.5) 
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In the 2x2 notation the above equation looks hke 

1(1 + 0^3) 



T 

n—l,n- 



-1(1 + 0^3) 







5(1 + ^3) 



1(1+^3); as) 
\ 







1(1 + ^3) 



(A.6) 



In the 4x4 notation the above equation becomes, 



T„_i,„-i = -[m(l + a,) + i5;„(7" - 7^') + PziYl' ' ll] , 



(A.7) 



From the matrix multiplication in the left hand side of Eq. ()A.1|) it can be seen that 



Tn-i,n is given as, 

^0 IN^ 


0-100 

Vo oy 

In the 2x2 notation the above equation looks like. 



Tn-i,n = ^2neB 



{k.i 



Tn-i,n = V2neB 







(A.9) 



^ + ia2j J 

Here ai and a2 are the first two Pauli matrices. When converted back to the 4x4 notation 
the above equation becomes. 



T„_i,„ = -^V2nei3(7i + 272) 



(A.IO) 



Similarly T„_„_i is given by. 



In the 2x2 notation the above equation looks like 



f ON^ 

10 



V-1 oy 



(A.ii; 



T„,„_i = V2neB I 1/ ■ n n 
\ -2(0-1 - ^0-2) 

which when converted back to the 4x4 notation becomes, 



cTi - 102) 



Tn-i,n = --V2nei3(7i - ^72) . 



(A.12) 



(A.13) 



Supplying the values of T^jS from Eq. (rO|) . Eq. (IXTfl . Eq. (|Xl0|l and Eq. ^KJ^ to 
Eq. ()A.1|) we get the result given in Eq. (j2.27|) . A similar procedure is followed to obtain 

Eq. ^rm . 
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APPENDIX B 



Phase factor of the Schwinger 

propagator 



B.l Loops of charged particles 

As we have discussed previously in chapter 3, the phase factor in the Schwinger prop- 
agator appears because the fermion propagator attaches two points with different gauge 
transformation properties. The Schwinger mechanism is best suited for processes where 
we do not have external charged particles and all the charged particles are inside the loop. 
For example one of the figures that contribute to the neutrino self-energy and another 
one which contributes for two neutrino three photon scattering as shown in Fig. 11.11 are 
perfect examples for these kind of diagrams. One important fact that can be noticed from 
the form of the propagators of charged fermions or gauge bosons [64j is that the form of 
the phase factor^^ of the different propagators are the same as given in Eq. ()3.1()|1 . To 
understand its importance we take concrete examples. 

First we take a one loop two point function containing virtual charged fermion internal 
lines. The one loop photon vacuum polarization diagram as shown in Fig. IB. II is a good 
example. The electron propagators connect points P and Q in space time. If we are 
interested in finding out the overall phase factor accompanying the vacuum polarization 
tensor then we will have to use the Schwinger propagator for the electrons. The contri- 
bution from the phase factors k{Q,P) and k{P,Q) which we denote as \1/(P, Q) will be 
according to Eq. (|3.6|) and Eq. (j3.1(J|) 



^^In the reference cited the authors did not explicitly take A(^) as they were not interested in the 
general properties of the phase factor. 



k{Q,P)k{P,Q) 
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Figure B.l: One loop two point function containing virtual charged fermion internal lines, 
Vacuum polarization diagram of photon. 

which reduces to unity because of the antisymmetry of F^j,. From Eq. ()B.1|) it is seen 
that the phase factor's contribution in the one loop calculation is trivial and obviously 
gauge invariant. In any kind of a loop where we have two charged particles the phase 
factor contribution is unimportant as all the charged particles have the same form of 
phase factor with the propagators. 

Next we take a loop with three charged particle propagators connecting three space- 
time points P, Q, R. We have shown that the phase factor between two points do not 
depend upon the path which joins them. Consequently to find the contribution of the 
phase factors in the calculation for the loop we join the three points by straight lines as 
shown in Fig. IB. 21 The overall phase contribution can then be calculated using Eq. ()3.6|) 
and Eq. ()3.10p and is given by 



As all the phase factors are of the same form the first point to notice is that the contribu- 
tion from the function A(^) cancel out in the overall factor, showing that the contribution 
is explicitly gauge invariant. 

The next point which requires to be discussed is about the path independence of the 
phase factor. From subsection 13.2.21 we know that k{Q, P) is independent of the path 
which joins them, but here the path independence of k{Q, P) does not imply 



^{P,Q,R) 



h{Q,P)k{R,Q)k{P,R) 




(B.2) 



^{P,Q,R) 



k{Q,P)k{R,Q)k{P,R) = 1, 



(B.3) 
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or from Eq. ()3.10|) 

/(P, P) = liQ, P) + liR, Q) + /(P, R) = 0. (B.4) 

Instead of the above expectation we get a finite contribution from Eq. ()B.2|) . As I{Q, P) 
consists of the product of the two end points instead of their difference, the different 
phase factors from the different paths connecting two intermediate points of the loop 
when multiphed does not reduce to unity. 

It can be shown in general that if we include loops where we have more than three 
charged particle propagators joining more than three points in space time then also the 
phase factor is going to contribute. To understand where we can actually neglect the 
phase factor and where not we have to calculate the flux of the uniform magnetic field 
attached to the area of the loop. 

B.2 The flux rule 

The phase factor contribution to the one loop calculations of various diagrams as 
shown in Fig. IB.l| and Fig. IB. 2l show that they are explicitly gauge invariant. In the case 
of Fig. IB. H it was shown that it contributes nothing for the phase factor. The contribution 
from Fig. IB. 21 can be understood in another way also. As the magnetic field is directed 
along the 2;-direction the component of area enclosed by the looping particles along the 
X — y plane is only active for its flux. From the path independence of k, we know that 
\I^(P, Q, R) can always be calculated from the triangular path connecting the three points 




Figure B.2: A loop consisting of three space-time points joined by charged particle propaga- 
tors. The arrow heads indicate the direction of charge flow. 
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Figure B.3: A loop consisting of four space-time points joined by charged particle propagators. 
The arrow heads indicate the direction of charge flow. 



as shown in Fig. lB.2l and so the flux attached with this loop must be the magnitude of the 
magnetic field B, times the area of this triangle in the x — y plane denoted by A±{P, Q, R). 
Calculating the flux we get 

^(P,g,i?) = BA^iP,Q,R) 

= ^-B[{Q-P)x{R-Q)l, (B.5) 

where [M x N]z represents z-component of the antisymmetric cross-product between two 
arbitrary 3- vectors M and N. Eq. ()B.5|1 can also be written as 

^(P, Q,R) = ^{Q- PfF^AR - QY , (B.6) 

and in this form we see from Eq. ()B.2j) that it can be written as 

^{P, Q, R) = exp {leJ^iP, Q,R)} . (B.7) 

Next we consider a loop containing four charged particle propagators. The four photon 
interaction in QED is a process where this loop occurs naturally. Here the the internal 
virtual particles may be electrons. In coordinate space the four points where the photons 
attach to the loop may be taken as P, Q, R and S. As far as the contribution from the 
overall phase is concerned we can join them by straight lines forming a general quadri- 
lateral as shown in Fig. IB. 31 Using the same technique for finding the overall phase as 
discussed in section IB. II we find here 

^(P,Q,i?,5) = k{Q,P)k{R,Q)k{S,R)k{P,S) 

= exp l^e^ [P^F^^Q'' + Q^F^^R" + R^F^.S" + S^F^.P''] | , (B.8) 
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which can also be written as 

^(P, Q, R, S) = exp {leJ^iP, Q, R, S)} , (B.9) 

where 

^(P, Q, R,S) = ^ HQ - PYF,,{R - QY + {S- RYF^,{P - Sf] . (B.IO) 

Working similarly as was done in the case of Eq. ()B.6|) it can be shown that the above 
flux is just the sum of the areas of the triangles PQR and SRP along the x — y plane 
times the magnitude of the magnetic fleld. 

From Eq. ()B.7|) and Eq. ()B.9|) we can generalize that the overall phase depends on the 
flux of the magnetic fleld attached with the area of the loop. Given an arbitrary diagram 
we can initially calculate this flux to flnd out the overall phase. Also this connection 
with magnetic flux explains another important fact that this phase is absent for a loop 
consisting of just two points. If we join this two points with straight lines then in no way 
will those straight lines enclose any area and as a result no resultant flux will emerge. 

The Schwinger propagator as given in Eq. ()3.2p is speciflc for the case of those gauges 
which gives a background magnetic fleld. In the general derivation by Schwinger the 
form of the fermionic propagator was derived for general constant gauge flelds. If we 
choose a pure gauge fleld which does not give rise to any electric or magnetic fleld as 

^^(0 = i^5/ea, (B.ll) 

where i^' is a constant, then also we can flnd the Schwinger propagator and it differs from 
the vacuum propagator only by the phase factor. Denoting this propagator in presence 
of a pure gauge as given in Eq. ()B.11|) as iSA{x, x') we have 

iSa[x, X = exp{tel{x, x / 7^7^ ^ , B.12 

where I{x,x') is given by Eq. ()3.13|) . From generalized Stokes theorem, which is 

fdeA^{0 = J da^'^F^^iO, (B.13) 

where da'^'^ is the inflnitesimal surface area in the — C,u plane and F^^, is the fleld 
strength tensor, we can say that I{x,x') is path independent as here P^,y = 0. The path 
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independence of /(x, x') allows us to join the two points by a straight line and proceeding 
similarly as done in chapter El while discussing about the phase factor, we get 

Jx' 

= K{x''x^-x"'x'J. (B.14) 

The above result implies that the overall phase factor for diagrams like those in Fig. IB. 21 
and Fig. IB. 3| will always be unity. As a result rates of those processes will remain similar 
to the vacuum results. 
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APPENDIX C 



A short note on 'Real Time 
Formulation' 

C.l A brief introduction to thermal propagators 

There are more than one way in which we can proceed in real time formahsm, out 
of them here we take the canonical approach [SEl as is prevalent in the zero temperature 
case. The main purpose of this section is to give a glimpse of the real time formalism, 
to show how a genuinely temperature dependent term and a temperature independent 
term appears in the propagators. Also some of the relevant formulas will be cited in this 
section. 

In this section, the essence of the real time formulation is presented in a heuristic 
manner using the real scalar fields. Suppose that we have a Lagrangian, 

^ = -d^cpd^<P-—<P\ (C.l) 

which describes a real scalar field. The Fourier mode expansion of the free field is given 

as, 

<P{x) = I , "^'^ VW^"" + a\vy'"^\ , (C.2) 

where the a{p) and a){p) are the annihilation and the creation operators of the bosons 
with momentum p. Formally the propagator is represented as 

^A^(x-i/) = (O|T0(a;)0(i/)|O), (C.3) 

where T is the time ordered product defined as 



T(t){x)(t){y) = Q{xo - yo)(t){x)(t){y) + 6(^0 - xo)(t){y)(t>{x) 
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(C.4) 



For obtaining the expression of the propagator from Eq. ()C.3|) we use the fact that 
(o I I 0^ = 0, always due to the basic definition of the operators a{p) and a"''(p), 

and get 

^^(^ -y)= — ■ (c-5) 

J [In) — 171'^ + le 

In a medium there are real particles and so the expectation value of the number operator 
a) a is not zero (a) {p)a{p)^ ^ ^ 0. So now we take 

a\p)a{p'))=5\p-p')fB{p), (C.6) 
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where 

1 



fB{p) 



is the Bose- Einstein distribution function. Here u is the four-velocity of the centre of mass 
of the heat bath. By the commutation relations naturally we have 

aip)a\p'))^ = 5%p - p') [fsip) + 1] . (C.7) 

In a thermal medium the following relations 

'^o) {p)o) {p')) ^ = and {a{p)a{p))p = 0, 

remain the same as those in vacuum. Utilizing these facts and doing the calculation as is 
done in the zero temperature case we find the expression of the thermal propagator as 

^f{p) = ^-Y—- - 2vr^5(p^ - m^)^-^ , (C.8) 

which has a genuine temperature dependent term as the second term in the right hand 
side. Here the temperature dependent term comes in with a delta-function showing this 
propagator corresponds to exchange of real particles unlike the exchange of virtual par- 
ticles present in the first term. The expression becomes realistic because the thermal 
bath always has some real particles, so in any process these particles can be exchanged. 
Eq. ()C8|) is really not the whole story, there are other propagators also accompanying 
it. In the next section some technical details is given which will elucidate the complex 
structure of this formalism. 
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C.2 A bit of formal theory 

C.2.1 The n-point Green's function in vacuum 

The Fourier expansion of the scalar field is given in Eq. ()C.2|) . When we are talk- 
ing about interacting field theory the free field expansion as given in Eq. ()C.2|) is often 
called To do interacting field theory we can always find an operator U which 

transforms the incoming free fields (pin to the interacting Heisenberg fields and the 

transformation is given by 



X) =U (j)in{x) U, 



t = x". 



For the above transformation to hold, we must have 



U{t) = Texp 



I I I Lb -^^intYViny 

— oo 



X] 



which is the evolution operator. It must satisfy the following property, 

lim U(t) = 1. 

Generally an n-point Green's function in the interacting theory is defined as 

G{xi, •, •, •,x„) = (O|T{0(xi)0(x2) ■ ■ ■ 0(x„)}|O). 
Introducing a very large time r such that 



5 5 5 5 *'n 



tn < r 



(C.9) 



(C.IO) 



(Gil) 



(G12) 



(G13) 



and choosing the convention 



U{T,-r)=Uir)U-\-T), (G14) 
we can manipulate the expression inside the time ordering. Using Eq. ()C.9|) . 



T{0(xi)0(x2) • ■ • 0(X„)} 

= T{f/-i(ti)</),„(a;i) U{h)U-\t2)<l)Ux2)U{t2) ■ ■ ■ U-\tn)<P{xn)U{tn)} 

= U-\t) [TU{T)U-\h)<P,n{Xi)U{h) ■ ■ ■ U-\tn)<P,n{Xn)U{tn)U-\-T)\ U{-r) 

= U'\t) [t{0,„(xi)0,„(x2) ■ ■ ■0,„(x„,)}t/-i(r,-r)] U{-r). (G15) 
Using the above equation we can write 

^(xi, x„) = (0|f/'i(r) [T{0,„(xi) ■ ■ ■ Mxn) U{t, -r)}] U{-r)\0). (G16) 
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In vacuum 



lim U(t)\0) = e'^lO) , 
lim f/(-r)|0) = 10) , 

T— ►OO 

where S is some real number. To take account of the phase the Green's function in vacuum 
is defined as 

^ (^i'-'-'-'^'^)- (0|f/(oo)|0) ' ^^-^^^ 

where 

f/(oo) = f/(r, -r) , (C.18) 

when T tends to infinity. The perturbation expansion is obtained by expanding f/(oo) 
and reducing the products of field operators through Wick's theorem. 

C.2.2 Interacting Scalar Fields in Presence of a Thermal Bath. 

When there is a thermal bath present, then 

lim f/(-r)|0,/3) = \0,P), 
(where |0,/3) is the thermal vacuum)but 

hm f/(r)|0,/5) ^e^^lO,/?), 

T — >00 

as there are real particles inside the thermal bath (thermal vacuum), they can interact 
and the state may change physically. As a result (0, /3|t/^^(r) cannot be just replaced by 
e-*'^(0,/3| as is done in Eq. (inTfll . Here 

G{xr, -, -, -, Xn) = {U-\oo) [T{<f)Uxi) ■ ■ ■ Mxn) U{oo)}]). (C.19) 
The expectation value of some operator O is defined as 

(O)^^. (CO, 

where Tr indicates a trace has to be taken (with respect to any basis) and Z is the 
partition function given by 



Z = exp 



-f3H + J2aAQA 

A 

90 



(C.21) 



where Qa are the conserved charges that commute with the Lagrangian and ua are 
the chemical potentials that parameterize the composition of the medium. Here we are 
working in presence of a thermal bath, and it can have a velocity. For the next part of 
the discussions we will work in the rest frame of the heat bath. 

The perturbation expansion is obtained by expanding the formal formula for U{oo), 
but in this case we have to take into account the factor U~^{oo), which is given by 



(C.22) 



U^^{oo) = Texp 
and expand it also. T indicates anti-time ordered product. 

C.2.3 Calculation of the Thermal Green's Function. 

Writing ^int[(pin{x)] as ^int{x) for brevity, the whole expression of the Green's function 
can be written as, 



G{X,, ; ; ; = £ 11 ^ {T \ f d'z, [ d% ■ ■ f d\Xnt{zi)Xntiz2) ■ -^Uz. 



X T 



(t>in{xi) ■ ■(t>in{Xn) j d^ ■ ■d\^int{yi) ■ -^intiyp) 



(C.23) 



Concentrating on (f)^ theory where the interaction term is given by 



we can write the Green's function as 

G{Xi, •, •, •, Xn) = 2^ j- 2^ : ( T 



s=0 

X T 



s\ 



p=0 



d'^zi I d'^Z2 ■ ■ / d'^z, 



4! 



4! 





s,p 


=0 




T{ 



n{Xl)- ■(t)in{Xn) j d^lji ■ -d^ 

i\y{-i\)P 



d z\ 

s\p: 



4! 4! 
d'^Zs / d'^yi ■ ■ ■ j d'^yp 



4! 



4! 



}T{(l)in{xi) ■ ■ ■ (piniXn)' 



4! 



^tniyp\y 



4! 



(C.24) 



This expansion shows that the general expectation value term will be of the form 

{T(AB ■ ■ ■F)T{AB ■ ■ ■ F)) (C.25) 
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where the bars over A,B,... indicate they are operators under the T sign, or they are 
anti-time ordered operators. The objects A,B,.. are the fields that appear in the last 
line of the Eq. ()C.24|) . The vertices that come under T are called the type-1 vertices, and 
those that come under T are called the type-2 vertices. A propagator which joins two 
vertices can now be of four kinds, 

• 1-1 type, which joins two vertices under the T sign. 

• 2-2 type, which joins two vertices under the T sign. 

• 1-2 type, which joins a type-1 vertex to a type-2 vertex. 

• 2-1 type, which joins a type-2 vertex to a type-1 vertex. 

Another important thing to be noticed is, any external vertex where real particles come 
in or go out must be of type-1 variety, as those fields 0(x)s are under the T operation in 
Eq. ()C.24|) . If now we calculate the two point function where both the fields are of type-1 
by simplifying the Eq. ()C.25|) properly then we will arrive at Eq. ()C.8|) . But there can be 
other options as well, where the fields are both of type-2, or are mixed. Not going into 
further details here the other propagators are supphed for the scalar theory. Apii{k) is 
given by Eq. ()C.8|) . the others are 



Af22(p) 



-1 



2TTiS{p^ — m'^)r](p ■ u) , 

w?')[ri{p ■ u) + 6{—p ■ u)] , 
m^) [rjiyp ■ u) + 0{p ■ u)] , 



(C.26) 



p- 



|2 — + ie 



Af2i(p) 
Afi2(p) 



2'Ki5{p^ - 
27ri(5(p^ - 



(C.27) 
(C.28) 



where 



rj{p ■ u) 



6{p ■ u)nB{x) + 6{—p ■ u)nB{—x) , 



(C.29) 



with 



1 



(C.30) 



e^' - 1 



and X is defined as. 



X = (3p ■ u . 



(C.31) 
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For fermions the procedure is same, except here we have anti-commutation relations 
in place of commutation relations, and the expressions for the four propagators are, 

1 



Spiiip) = (^ + m) 



Sf22{p) ^ (J) + m) 



— + 2'KiS{p -m)rjF{p-u) 
— m-' + le 

-1 



„ „ .- + 2m5{Tp' - rn^)r]Fip ■ u) 
— + i£ 

SfMp) = + m)2TTi6{p'^ - ■m'^)[r]F{p ■ u) - 9{-p ■ u)] , 
Sf2i{p) = {.i> + m)2TTi5{p^ — m'^)[r]F{p ■ u) — 9{p ■ u)] . 

In the case of fermions r]F{p • u) is given by, 

r]F{p ■ u) = e(p • u)fF(p, (3) + 0(-p ■ u)fF(-p, -//, (3) , 



where. 



fF{p, 1^,(3) 



(C.32) 

(C.33) 

(C.34) 
(C.35) 



(C.36) 



(C.37) 



e^(p-u-/i) _|_ I ' 

is the Fermi-Dirac distribution function. 

With this brief exposition of statistical field theory in the real time formalism we 
understand why instead of one propagator we have four of them. But usually most of 
the time we only work with the 1-1 component as is done in this thesis. But formal 
calculations in statistical field theory requires all four propagators as specified above. 
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APPENDIX D 



Manipulations under the integral 

sign 



All the results appearing in this appendix have been derived previously in (SUl- They 
are explicitly written down here because the following techniques are very useful in ma- 
nipulating the integrand inside the integral sign in an equation like Eq. ()4.59j) . If we 
concentrate on the rest frame of the medium, then p ■ u = Po- Thus, the distribution 
function does not depend on the spatial components of p. From the form of Eq. (j4.6(J|) 
and Eq. (|4.61|) we find that, the integral over the transverse components of p has the 
following generic structure: 

J (fp^ e*(^''^)e*(P''^') X {pP^ or p'^^) . (D.l) 

Notice now that 



dp0 



^ ' ^ftan(ei3s)/^ + tan(eSs')p'^^)e*(^''^)e*(P''^'). (D.2) 



g*(p.«)g*(p',s') 



eB 

However, this expression, being a total derivative, should integrate to zero. Thus we 
obtain 

J d^p^ tan(eSs)/^e*(P'^)e*(^'''^') = ~J d^p^ tan(eSs')p'''^e*(^''^)e*(^'''^') . (D.3) 

This above equation is symbolically written as, 

tan(ei3s) = - tan(e;Bs') p'^^ , (D.4) 

where the sign '=' means that the expressions on both sides of it, though not necessarily 
equal algebraically, yield the same integral. This gives 

. ^ _ tanjeBs') 

^ tan(e;Bs) + tan(e;Bs') ' ^ ' ' 

,n o tan(ei3s) , ^ s 

^ tan(ei3s) + tan(ei3s') ' ^ ' ^ 
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Similarly we can derive some other relations which can be used under the momentum 
integral signs. To write them in a useful form, we turn to Eq. ()D.2|) and take another 
derivative with respect to •p'^^ . From the fact that this derivative should also vanish on p 
integration, we find 

1 



a /3 _2_ 

^^■^^ ~ tan(ei3s) + tan(ei3s') 
In particular, then, 

1 



ieS af3 tan^(ei3s') 

^ tan(ei3s) + tan(ei3s') 



1% \ Hj \ 



(D.7) 



2 o 
P± = 



tan(ei3s) + tan(ei3s') 
It then simply follows that 

p± - 



-ieB + 



tan(ei3s) + tan(eSs') 



-ieB + 



tan(ei3s) + tan(ei3s') 



tan^(ei3s) 
tan(ei3s) + tan(ei3s') 



(D.8) 



(D.9) 



And finally using the definition of the exponential factor in Eq. ()H.4|1 we can write 



m 



2 o 



+ (pfi — sec^(e;Bs)p^ 
. as " 



(D.IO) 



Here the '=' means that both sides of the above equation, when multiplied by e'^'-^'^^ e'^^P' 
are equivalent inside the momentum integrals. 

This set of relations help us considerably to handle the terms present inside the inte- 
grals appearing in Eq. ()4.6H) and Eq. ()4.60|) . 
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APPENDIX E 



A list of relevant integrals 



These kinds of integrals often crop up when we are calculating with the Schwinger 
propagator. 



Here a is a real number. Obviously the second integral can be obtained from complex 
conjugating the first one. 

Next we consider integrals containing vector indices. 



Here, as before a is a real number, p and B are 4-vectors. For each one of the above 
integrals there exists one corresponding complex conjugate integral. The set of complex 
conjugated integrals is not written here. 

Finally comes those integrals which involve the parallel and perpendicular components 
of the vectors. 




(E.l) 



(E.2) 




(E.3) 



(E.4) 



(E.5) 




(E.6) 



(E.7) 



(E.8) 



(E.9) 



(E.IO) 
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Here also we have one corresponding complex conjugated integral for each one of the 
above integrals which are not written. 
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APPENDIX F 



Some calculations in a magnetized 

medium 



F.l Derivation of Rjfj and Rjfj 

In this appendix those calculations are done which were side stepped in chapter |1J 
First it is shown how we get Eq. ()4.65|1 from Eq. ()4.60|1 and how Eq. ()4.66j) is derived. 
From Eq. ()4.60|1 it is clear that if we apply Eq. (jP.Tjl and then Eq. ()D.5|) to the last term 
in the right hand side then it drops out and we get Eq. ()4.65|1 . 

An expression of R^°J is given in Eq. ()4.61|) . leaving it aside we derive Eq. ()4.66p in a 
different way. Now from Eq. (|4.59p it is known that in the rest frame of the medium, 

^t,uip,p\ s, s') = Tr[7^75G(p, s)-f^G{p, s')]??f(po) + Tr[7^75G(-p', s')-f^G{-p, s)]T]Fi-po) ■ 

(F.l) 

R^!y(p, p', s, s') contains the sum of two traces multiplied by the functions riF{±po). Except 
the 7] functions, the two traces are equal, so for the part which is odd in orders of the 
external field, Eq. ()F.1|) is equivalent to, 

R!rJ=^+MTr[T(^)] (F.2) 

where for brevity the arguments of R^°} and T^°^ are suppressed. Here rj+ipo) = i]{po) + 
vi~Po) and T^°^ is given by, 

T^l^ = ztan(ei3s) rn^-f^-f^azju + lirih<yzi)\\lvi>\\ + 7m750-2|^||7!^/± sec^(e;Bs') 
+ 2tan(ei3s') [m^7M757vO-2 + lia^i>\\"i^^^i>\\ + IvH^i-lvOz^w sec^(e-Bs) 

(F.3) 
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Writing out the terms in Tj°\ we get 

+ ip'e-p'k'h,^,] , (F.4) 

Tr [7^751^11 7i^cr^l^||] = Tr [2p,.||7^75cr^|^j| + pl'J^.'Julbf^z + {p ■ k)\i'j^'yul5(^z 

- (p°A:=^-p'fc°)7/.7.] ■ (F.5) 
Using Eq. ()F.4|) and Eq. ()F.5|) we can write Eq. ()F.3|) as, 

Tr[T^°)] = ztan(ei3s) Tr [{m' - pj - ■ k)^^}j,-f^-f,a, + {p^k' - p'k'>)-f,-f^ 
+ Vj^ii 7^750-^1^11 +7^75a;,|^||7i./_Lsec^(e;Bs') 

+ itan{eBs') Tr [-{m^ -Py - (p • A;)||}7^7,.75cx^ - (/A;^ -p^k^hf.lu 

+ 2pi,||7^75a-^|^jj + 7^75|^_l7i^o-2|^|| sec^(ei3s) . (F.6) 

As we are interested in taking the trace of Tjf), so now we write down the traces of the 
components of Tj^^. The easier ones comes first, and they are 



Tr[7^75a^|^j|] = 4:6 ^l2a^^p"'^^ , 
Tr[7^75a^j^l|] = 4e^i2a||P'"ii , 
Tr [7^7^,75(7^] = - Tr [7^7^75(7;, 
= —4e^yi2 ■ 

Next comes the trace of a part of T^J ^ and it is 
Tr ' 



(F.7) 
(F.8) 

(F.9) 



lpnh(yzi>\\lvi>^ tan(ei3s) sec {eEs) + ^pn^i)invOzi>\\ tan(ei3s') sec {eEs) 
= ^9fMa^\ sec^(ei3s') tan(ei3s) + p"^^^Puj^ sec^(ei3s) tan(ei3s') 

+ 45(^„|| sec^(e;Bs') tan(e;Bs) + sec^(eSs) tan(eSs') 

+ 4c/^„||fc°'iip^^ sec2(ei3s)tan(ei3s') . (F.IO) 

Now using Eq. ()D.5|) and Eq. ()D.6|1 from appendix [Ul the above trace can be rewritten 



as, 



Tr 



7fil5<^z^\\lu^'± tan(eSs) sec^(ei3s') + ^^j.^5^±^,ycrzfn tan(ei3s') sec^(ei3s) 



= ^{g^wP'^^^kyi^ +5',.Q|^p"iiA;^^){tan(ei3s) -tan(ei3s')} 
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Combining the previous steps, the trace of T^J can be written as, 

Tr[r^°)] = Ai [5^,12 {m^ - - (p ■ A;)||}{tan(ei3s) + tan(ei3s')} 
+ 2£^i2a||K,|P"ii tan(ei3s) tan(ei3s')} 
+ g^.u{p^k^ - P^A;°){tan(eSs) - tan(ei3s')} 
+ P'^Kdt^a^^K^ + gua\\k^^){iim{eBs) -tan(ei3s')} 

~ sec^(ei3s) tan^(ei3s') ~ 2/ v, \ 1 v-> i\ 

+ ^-"^ "^-^ tan(eig.)+tan(eig.O ^'^"" "'''^^ (e^^) tan(eS. ) 



(F.12) 



Now we modify the term [w? — p|){tan(e;Bs) + tan(ei3s')} in the first hne of the right 
hand side of the above equation. As the trace of T^^ is related to R|fj by Eq. ()F.2|) . 
and as R^°j is inside the momentum integrals in Eq. ()4.59|) in chapter |3 so we can use 
the operations under the integral signs. The operations were explicitly written down in 
appendix |D| Using Eq. ()D.10|) we can write. 



-J9j){tan(ei3s) + tan(eSs')}e[*(^''^)+*(f''^'W 
= {tan{eBs) +tan(ei3s')} 



i- sec^(eSs)p^ 

as 



,[$(p,^i)+*(p',s')] 



where e'**^^''^^^**^^'''*'^] is the overall phase factor appearing inside the momentum integrals 
in Eq. ()4.59p . Using Eq. ()D.8|1 the above equation becomes. 



— Py){tan(e;Bs) + tan(eSs')}e 



:<i>(p,s)+$(p',s')] 



{tan(ei3s) + tan(ei3s')}(i— 

as 



sec2(ei3s) <^ -ieB + 



tan2(ei3s') 



tan(ei3s) + tan(ei3s') 



J*{p,s)+$(p',s') 



(F.13) 



The above equation is valid under the momentum integrals and the parametric integrals in 
Eq. ()4.59|) . so we can do some of the parametric integrals right now to modify Eq. (jEH- 
First of all, 



tan(ei3s');f e[*(^'^)+*(p''^'» ds = tan(ei3s')e*(P''^') 

= 0, 



(F.14) 



then. 



tan(ei3s)-^e[*(P'^)+*(^''^')l ds = -leB / sec2(eSs)e[*(^'^)+*(P''^')l ds , (F.15) 
00 us J— 00 
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and so combining the above steps we get from Eq. ()F.13|) . 



p^){tan(eS.)+tan(ei3.')}et*(^'^)+*(^'^')l = _ ^^^'(f ^) ^^^'(^^^0 
^11^^ ^ ' ^ tan(ei3s) + tan(ei3s') ^ 



X g[*(p.'5)+*{p'.s')] _ (F.16) 

Now from Eq. (lFl2l) and Eq. (IRTgI) . we can write 
Tr[T(°)] = 4z 



f sec^ (ei3s) tan^(ei3s') ,9 , , , ^ x / ^ 

e^,ul2 { ( u + ■p)||{tan(eSs) +tan(ei3s')} 



I tan(ei3s) + tan(ei3s 
+ 2e^i2a||(p'j,||]9"" tan(ei3s) + Pu^^p"^^^ tan(ei3s')) 

i 7: r /^/N-, sec^(eSs) tan^(eSs') 

+ ^M"!!^'^^ "{tan(eSs) -tan(e;Bs')} - " v ; v ; 



tan(ei3s) + tan(ei3s') j 
+ {9,iu{p ■ k)\\ + fi'i.a||P""fc/.±}{tan(eSs) - tan(ei3s')} 

+ gua^^k^^^P^^±sec^{eBs)tan{eBs')]. (F.17) 

From the above equation and Eq. ()F.2|) we immediately see that R^°j matches with 
Eq. (pr66|l in chapter H 

F.2 Derivation of the expression of IIIq 

First of all it is shown why the other components of U^^q except IIIq turns out to be 
zero in the limit when the external momentum goes to zero. First the contributions from 
nj^o'^ are considered. From Eq. ()4.59|1 and Eq. ()4.fjfjj) we can write 



{2ny 

c/s'e*(P''^')r/+(po) [2eou3{py tan{eBs) + pop'Han{eBs')} 
+ ^oo(p- ^)||{tan(ei3s) -tan(ei3s')}l • (F.18) 



X 







jg g^j^ even function of po when = 0. ?7+(po) inside the integrand is also 
an even function of po- From these observations it can be said that when ko = the 
integrand on the right hand side of the above equation becomes an odd function of po and 
result the integral vanishes. 
Then comes H^q"'' and it is, 



(2vr)^ J-c 

X / ds'e*(^'''^')r/+(po) p^k^{tan{el3s) -tan{eBs')} + kYsec'^ieBs)tan{eBs') 

Jo L 

(F.19) 
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n^o can exist only when the parametric integrals can yield terms which cancels terms 
like or k^. In the present case as this does not happen so due to the presence of the 
terms like k^ and k^, the integral vanishes in the limit fc — > 0. Similarly 1120°^ will also be 
zero. 

Now comes the contributions from n^[f\ From the expression of in Eq. (lOKl) we 
can immediately say that TI^J^^^q must be zero. If on the other hand /i = 1, then 

X / ds'e*(P''^')r/+(po)£ioa/3{p""p''^^sec2(ei3s') +p°V^"sec2(ei3s)}, (F.20) 
Jo 

and the above equation shows that each term in the right hand side inside the momentum 
integrals contain one perpendicular component of loop momentum, and so it must inte- 
grate out to perpendicular components of the external momentum k. H^q*^^ can exist when 
the parametric integrals can produce terms which will cancel the perpendicular compo- 
nents of fc. As in the present case it does not happen as a result the integral vanishes 
when fc — s> 0. Similar result is obtained for 1^20^^ ■ Ultimately the only component of II^q 
which is non vanishing in the specified external momentum limit is Ilgg^^ which is simply 
written as II^q. 

Now the relevant steps are given by which we get Eq. (j4.75p from Eq. (j4.74p . To begin 
with we explicitly write down the sum of the functions s) and in terms of 

the momentums p and p'. 

s) + s') = isp\ + is'p'i^ - i{s + s')m'^ 

- ■^{iim{eBs)p\ + tan(e-Bs'y_^} - e\s\ - e\s\ . (F.21) 

Now we concentrate on the integration of the perpendicular components of the loop 
momentum in Eq. ()4.74j) . and try to see what happens when \k_^_\ ^ Q. It must be 
noted that we have thrown away some of the terms in deriving the expression of II^q by 
assuming that fco = and ^ 0. The terms still existing in Eq. ()4.74|) do not vanish in the 
above specified momentum limit. First we integrate out the perpendicular components of 
momenta in the integral appearing in Eq. ()4.74j) . As the integrand in the right hand side of 
Eq. ()4.74|) does not contain any term involving the perpendicular momentum components, 
except the phase like terms, so here the integration is straight forward. From Eq. ()F.21|) 
if we only take the terms containing the squares of the perpendicular components of p 
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and p', and try to integrate them out, using the results given in appendix 10 we will get 



i2n) 



2 exp { —[taii{eBs)p']_ + tan(ei3s')p'_f]} 



tan(ei3s) tan(ei3s') 
tan(ei3s) + tan(ei3s') 



} 



/ 777^ *5^P I ^[tan(ei3s) + tan(ei3s') 



P± + kj 



tan(ei3s') 



tan(e;Bs) + tan(e;Bs') 



(F.22) 



ieB 



and as k'j_ —>■ the above equation simplifies to, 

/ , ^."i exp I —\ta.n(eBs)p^< + ta.n(eBs')p'f]} = , ^ , — — , ^ ,, .(F.23) 

J (27r)2 eS^ ^ '^^ ^ '^^^^ 47r tan(ei3s) + tan(ei3s') ^ ' 

Using the above equation we get Eq. ()4.75j) from Eq. ()4.74|1 . 



F.3 Effective charge calculation 

This section gives the details of the calculations by which the expression of neutrino 
effective charge was calculated in chapter El 
The s integral in Eq. (j4.75p gives 



ds e 



is{pf,—rrP)—£\s\ 



(F.24) 



and the s' integral gives 



ds e 



/ is'(p''^—m^)—s\s'\ 



(F.25) 



— m^) + ie 

Using the above results in Eq. ()4.75|l and using the delta function constraint, we arrive at. 



'^pI 



(fc|f + 2(M)!|) 2 



(F.26) 



In deriving Eq. ()F.26jl . pieces proportional to A;j| in the numerator were neglected. Now if 
one makes the substitution, p',, [pu + ku/2) and sets k^ = one arrives at. 



+ 



P^h 2E' 



(F.27) 



Here n±{E'p) are the functions fpi^E'^, [3), and fF{Ep,—fi,P), as given in Eq. ()3.19p . 
which are nothing but the Fermi-Dirac distribution functions of the particles and the 
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antiparticles in the medium with a modified energy E'^. The new term E'^ is defined as 
follows, 

E;' = [ips-h/2)r + m^ 
and it can be expanded for small external momenta in the following way 

Ep c^pI + w? - p3k3 = El- psh 



where E^ = pl + rn? . Noting, that 



E'=Ep-^^ + 0{J4), (F.28) 



2Ep 



one can use this expansion in Eq. ()F.27|) . to arrive at: 

dpz 



= lim 2ie'B) I (n+(E')+n(E'- 



Ep 



In the classical gas limit the expression for for ri^{E'^) = n+{E'p) + n_{E'p) when expanded 
in powers of the external momentum is given by 



vAEl) = {l + l^)v4Ep), (F.30) 



up to first order terms in the external momentum k^. 
From Eq. and Eq. we get, 

dps ( Ep ^ (3' 
(27r)2 {p^k; 

The first term in the integral on the right hand side of the above equation contains a 
factor of p3 in the denominator and so is an odd function of p^. As a result the integral 
vanishes by symmetry. But the second integral survives and IIIq is given by, 

5 e^B ^ f dp 



nio = ^/5 J ^V+iEp). (F.32) 



In the classical limit: 

n^o = 0/?cosh(/?/i) J dp e-^^^ (F.33) 

Using the standard result 



g-aVb^+x^ cos(cx) dx = Ki{bVa^ + c^) , (F.34) 

Va^ + 

where a, b, c are real constants here, and putting c = 0, IIIq becomes, 

= —mp cosh{p fi)Ki{mp) . (F.35) 



TT 



Ki{m(3) is the modified Bessel function. 
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APPENDIX G 



Gauge invariance 



The axialvector-vector amplitude has an electromagnetic vertex and as a result the 
electromagnetic current must be conserved in that vertex. In our case the //-vertex is 
the electromagnetic vertex and the current conservation condition is given in Eq. ()4.Hfjj) . 
The following sections explicitly verifies the current conservation condition, which is also 
called the gauge invariance condition here. 

G.l Gauge invariance for 11^^, to even orders in the external field 



The axialvector-vector amplitude even in the external field is given by 

ililf = -e'J ^ jyse''^^-^ f ds'e''^^'-'^R^;}{p,p',s,s'). (G.l) 

Noting that, 

we can write Eq. ()4.65p as, 

R^] = 4^r/_(po) [{e^^appy^ - e^uap^py^^ - e^ua^pp''^p'^){l + tan(ei3s) tan(eSs')) 
+ £M-a/3xPV^^ sec\eBs') + sec\eBs)\ . (G.2) 

Here throughout we have omitted terms such as e^va^p^p'^^p'^^ i since by the application 
of Eq. ()D.5|1 we have 

'^lJ.ua±f3±P P — '^fiva^P^P P \ '^fj.ua^P^P 1^ 

SL tanjeBs') ^^^^ 

tan(eSs')+tan(eSs')^'""^" ' 

which is zero. 
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After rearranging the terms appearing in Eq. ()G.2|) . and by the apphcation of Eq. ()D.5|) 
and Eq. ()D.6|) we arrive at the expression 



B tiuai3P°' + tan(ei3s) tan(ei3s')) 

(G.3) 



....... .n.^.ta„(eB.).a„(eB/) ;-;:^;;;;-;g: 

Because of the presence of terms hke e^yapk^ and e^yai3±_k'^ if we contract R^fj by k^ , it 
vanishes. 

G.2 Gauge invariance for 11^^^ to odd orders in the external field 



The axialvector-vector amphtude odd in the external field is given by 

, ,74, 

9 1'- 



nj(;) = -e^ / ^ ds e*(^'^) f ds' e*(^''^')R(°) (p, p', (G.4) 

where R^°J(p, p', s, s') is given by Eq. (j4.66p . The general gauge invariance condition in 
this case 

fc'^nJ(;) = o, (G.5) 

can always be written down in terms of the following two equations, 

rnj;°) = , (G.6) 

rnjj) = , (G.7) 

where n^*j°j is that part of H^^^f^ where the index /i can take the values and 3 only. 
Similarly Il^J;fl stands for the part of H^^f^ where fx can take the values 1 and 2 only. H^j^j 
contains R|j°^y(p, p', s, s') which from Eq. (j4.66|) is as follows, 



Rj;), = 4zr^+(p 

+ 2e^||i2a|| tan(e;Bs) +Pv||P'"" tan(ei3s')) 

sec^(ei3s) tan^(ei3s' 



+ 9,,.,k^. (tan(ei3s) - tan(eSs')) - k^ ^^^^^^^^ ^ ^ 

+ gf^^^u{p-k)\i{tan{e8s)-tan{el3s'))] , (G.8) 

and n^^°^ contains R^°|^(p, p', s, s') which is 

R-^rli/ = 4i77+(po) {gtij_iy{p ■ k)\\ + 5f;,a||P"iiA;^^}(tan(ei3s) - tan(ei3s')) 

+ ^^.^iifc^iip^^ sec^(e-Bs)tan(eSs')] • (G.9) 
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Eq. ()G.6|) and Eq. ()G.7|) implies one should have the following relations satisfied, 



(G.IO) 



and 



-oo ' 



(27r)4 7-00 JO 

Out of the two above equations, Eq. ()G.10|) can be verified easily since 



(G.ll) 



k ^^^u — 0. 



(G.12) 



Now we look at Eq. ()G.11|) . We explicitly consider the case /xy = 3 (the /xy = case 
lead to similar result). For /iy = 3 



FRg = -Po 



{p'w ~ Pjj)(tan(ei3s) + tan(ei3s')) — kj_(ta.n{eBs) — tan(ei3s')) (4i?7+(po))- 

(G.13) 



Apart from the small convergence factors. 



^ s) + s')) = {p\f +pI- 2m') ^ - {p\f - pj) C - P'l tan(e - C) 



- pI tan(e + C) 



(G.14) 



where we have defined the parameters 

e = ^eBis + s') , 
C = leB{s - s') . 

From the last two equations we can write 

^eB ±e^M+Hp',s') ^ ^Hp,s)+Hp',s') ^^^2 _ ^2 _ ^^^2^^ - Q + pi sec\i + C) 
etc, 



(G.15) 



(G.16) 



which implies 



p'\\ ~ P\\= ^^^^ ~^ sec^(eSs') — sec^(ei3s) 



(G.17) 
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The equation above is valid in the sense that both sides of it actually acts upon e*'-^'*'^'''^'^ 
inside the momentum integrals, where 



s, s') = ^{p, s) + s' 



(G.18) 



From Eq. ^TJ^ and Eq. (|m7|l we have 



k^Rs^ = -4ir/+(po)po {p'^ sec^eBs') - pi sec^(ei3s))(tan(ei3s) +tan(e-Bs')) 



- A;i(tan(e;Bs) - tan(eSs')) + ieB{tan{eBs) + tan{eBs')) 



d_ 



(G.19) 



Now using the the expressions for pj_ and p'f from Eq. ()D.8|1 and Eq. ()D.9jl we can write 



r Rg^ = 4eBr]+{po)po [{sec\eBs) - sec\eBs')) 
+ (tan(ei3s) +tan(eSs'))^ 



(G.20) 



Explicitly writing e* in its proper place, the above equation can also be written as 



FRg^e* = 4:eBr]+{po)po-^ 

dQ 



e*(tan(ei3s) +tan(ei3s')) 



(G.21) 



Transforming to ^, Q variables and using the above equation we can write the parametric 
integrations (integrations over s and s') on the left hand side of Eq. ()G.11|) as 



oo roo 



ds / (is'^Rg^e" 

oo Jo 



eB 



oo POO (j 

dd rfCe(e-C);T7^(e,C) (G.22) 

-oo J —oo (2C, 



where 



Q = e*(tan(ei3s) + tan(ei3s')) • 



(G.23) 



The integration over the ^ and ( variables in Eq. ()G.22|) can be represented as. 



/oo /"OO rj 

dU rfce(e-c)-:^(e,C) 
-oo J~oo CLL 



-oo 
/•oo 



-oo 
roo 



/oo roo 
dC / dC 
-oo J —oo 

d^H^,0 



dC 



(G.24) 
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here the second step follows from the first one as the first integrand containing the G 
function vanishes at both limits of the integration. The remaining integral is now only a 
function of ^ and is even in pq. But in Eq. ()G.22|) we have ?7+(po)Po sitting, which makes 
the the integrand odd under po integration in the left hand side of Eq. (jd.llll . as f]+{po) 
is an even function in pq. So the po integral as it occurs in the left hand side of Eq. (jd.lljl 
vanishes as expected, yielding the required result shown in Eq. ()G.6|) . 
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